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Abstract 

We conclude the study of the Post-Minkowskian linearization of ADM tetrad gravity in the 
York canonical basis for asymptotically Minkowskian space-times in the family of non-harmonic 3- 
orthogonal gauges parametrized by the York time 3 K(t, a) (the inertial gauge variable, not existing 
in Newton gravity, describing the general relativistic remnant of the freedom in clock synchroniza- 
tion in the definition of the instantaneous 3-spaces). As matter we consider only N scalar point 
particles with a Grassmann regularization of the self-energies and with a ultraviolet cutoff making 
possible the PM linearization and the evaluation of the PM solution for the gravitational field. 

We study in detail all the properties of these PM space-times emphasizing their dependence on 
the gauge variable 3 /C(i) = ^ 3 -^(i) (the non-local York time): Riemann and Weyl tensors, 3-spaces, 
time-like and null geodesies, red-shift and luminosity distance. Then we study the Post-Newtonian 
(PN) expansion of the PM equations of motion of the particles. We find that in the two-body 
case at the 0.5PN order there is a damping (or anti-damping) term depending only on 3 /C(i)- This 
open the possibility to explain dark matter in Einstein theory as a relativistic inertial effect: the 
determination of 3 /C(i) from the masses and rotation curves of galaxies would give information 
on how to find a PM extension of the existing PN Celestial frame (ICRS) used as observational 
convention in the 4-dimensional description of stars and galaxies. Dark matter would describe 
the difference between the inertial and gravitational masses seen in the non-Euclidean 3-spaces, 
without a violation of their equality in the 4-dimensional space-time as required by the equivalence 
principle. 
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I. INTRODUCTION 



In Refs.fl, 2], quoted as papers I and II respectively, we studied Hamiltonian ADM 
tetrad gravity in asymptotically Minkowskian space-times in the York canonical basis defined 
in Ref.[3] and its Hamiltonian Post-Minkoskian (HPM) linearization in a family of non- 
harmonic 3-orthogonal gauges. Since in this formulation the instantaneous 3-spaces are well 
defined, we have control on the general relativistic remnant of the gauge freedom in clock 
synchronization, whose relevance for gravitational physics will be investigated in this paper, 
where the matter consists only of N scalar point particles (without the transverse electro- 
magnetic field present in papers I and II), in the Post-Minkowskian (PM) approximation. 

The definition of 3-spaces, a pre-requisite for the formulation of the Cauchy problem 
for the field equations, is done by using radar 4-coordinates a A = (o~ T = r;a r ), A — r,r, 
adapted to the admissible 3+1 splitting of the space-time and centered on an arbitrary 
time-like observer x M (r) (origin of the 3-coordinates a r )\ they define a non-inertial frame 
centered on the observer, so that they are observer and frame- dependent. The time variable 
r is an arbitrary monotonically increasing function of the proper time given by the atomic 
clock carried by the observer. The instantaneous 3-spaces identified by this convention for 
clock synchronization are denoted E T . The transformation a A h-> x^ = z^(t,o~ t ) to world 
4-coordinates defines the embedding z^{r,a) of the Riemannian instantaneous 3-spaces E r 
into the space-time. By choosing world 4-coordinates centered on the time-like observer, 
whose world-line is the time axis, we have x M (r) = (x°(r);0): the condition x°(t) = const. 
is equivalent to r = const, and identifies the instantaneous 3-space E r . If the time-like 
observer coincides with an asymptotic inertial observer x M (r) = x^ + e^r with = (1; 0), 
e£ = (0; 5 l r ), x% = {x°\ 0), then the natural embedding describing the given 3+1 splitting of 
space-time is z^{r, a r ) = x^ + e A a A and the world 4- metric is 4 g^ u = e A a Qab (e^ are flat 
asymptotic cotetrads, e A e% = 5 A , e A e A = 6%). 

From now on we shall denote the curvilinear 3-coordinates a r with the notation a for the 
sake of simplicity. Usually the convention of sum over repeated indices is used, except when 
there are too many summations. 

The 4-metric 4 ^ab has signature e (H ) with e = ± (the particle physics, e = +, and 

general relativity, e = — , conventions). Flat indices (a), a — o,a, are raised and lowered 

by the flat Minkowski metric 4 ?7( a )(/3) = e(H ). We define 4 ^( a )(6) = — e ^(a)(b) with a 

positive-definite Euclidean 3-metric. On each instantaneous 3-space E T we have that the 
4-metric has a direction-independent limit to the flat Minkowski 4-metric (the asymptotic 
background) at spatial infinity A g AB (T,ff) ->• 4 i]AB(asym) = e(H )• 

After a review of the York canonical basis and of the HPM linearization in Subsections 
A and B respectively, we will outline the new results of this paper in Subsection C. 

A. The York Canonical Basis 

In the York canonical basis of ADM tetrad gravity of paper I 
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the family of non-harmonic 3-orthogonal gauges is the family of Schwinger time gauges where 
we have 



a(o) (r, a) « 0, <p (a ) (r, <?) « 0, 



^(r,a) «0, tt^t,*) = -^*K(t,#) « ^F(r,ff), (1.2) 

where -F(r, a) is an arbitrary numerical function parametrizing the residual gauge freedom in 
clock synchronization, namely in the choice of the non-dynamical aspect of the instantaneous 
3-spaces E T . 

oA o^4 q a 

In the York canonical basis we have ( 4 E^ are arbitrary tetrads; 4 E^ and 4 -E( Q ), 4 -^( a ) = 

-R(a)(fe)( a (e)) 4 -^(b)! are tetrads adapted to the 3-spaces; 3 e T ( a ), 3 ^( a )r an d 3 e( a - ) , 3 e( a )r, with 
3 e( a) = i?( a )( b )(a(e)) 3 e[ 6) , 3 e( a ) r = i2( a )( b )(a( e )) 3 e( 6 ) r , are triads and cotriads on the 3-spaces 
S T ; for the shift function we have ri( ) = i2( tt )(6)(a( e )) n( a )) 

A A 
V)(y?( C )) + 4 ^( 6 ) i^)(a)(«(c)) ^ (a) ( Q )(^(c)) « 4 ^(o), 

(o) = 4 ^( o) = ^ (i; -n (a) %)) = / A , 4 %) = (0; %)) - %), 

(l + n)(l;0) = e/ A , 4^ = (n (o) ; 3 e (a)r ) « 4^, 

o (a) o 08) o (a) o (/?) 

4_ _ 4t7i(Q!)4^ 4 771 (P) _ 4 771 4_ 4771 —4771 4_ 4t7> 
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g Tr = ^(a) 3 e(a)r = ~e(j) 1/3 Q r , 

a 

[ 9rs = -€ 3 g rs = -e 3 e(a)r 3 e( a ) s = -e0 43 ^r S = -e <fi 2/3 Q 2 r . 5 rs , 



a 

Q a = e r " 1) = e E * 2 t*-^, = 6 = ^7 = y/det^g = 3 e, 

3 e ( a)r = V3 5ra, 3 e[ a ) = V 1/3 <J ra . (1.3) 

The set of numerical parameters 7« a satisfies [3] Xl« 7a« = 0, J2 U 7a«75« = ^a5, 
7a«7af — — §■ A different York canonical basis is associated to each solution of 
these equations. Let us remember [4] that to avoid coordinate singularities we must always 
have N(t,<j) = 1 + n(r,a) > (3-spaces at different times do not intersect each other), 
e A g TT (r, a) > (no rotating disk pathology) and 3 g rs (r, <?) with three distinct positive eigen- 
values (M0ller conditions). 
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B. The HPM Linearization 

The standard decomposition used for the weak field approximation in the harmonic gauges 

is 

9[tv V^v h[tv 5 | h^i/ 1 , 1 a h^ v | , 1 a Op h^ v | < < 1 , 

(1.4) 

where ^n^ is the flat metric in an inertial frame of the background Minkowski space-time. 
This is equivalent to take a 3+1 splitting of our space-time with an inertial foliation, having 
Euclidean instantaneous 3-spaces, against the equivalence principle and against the fact 
(explicitly shown in paper I) that each solution of Einstein's equations has an associated 
dynamically selected preferred 3+1 splitting. The final x° = constant 3-spaces have a non- 
zero extrinsic curvature but are not well defined in this formalism and depend on the initial 
data at r — > — oo needed to solve the wave equations in the linearized harmonic gauges for 
the lapse and shift functions (see paper II). 

Instead the HPM-linearization of paper II of the Hamilton-Dirac equations in the (non- 
harmonic) 3-orthogonal Schwinger time gauges (1.2) uses as background the asymptotic 
Minkowski 4-metric existing in our asymptotically Minkowskian space-times. By using radar 
4-coordinates adapted to an admissible 3+1 splitting of space-time, we put 

4 g AB (T,<j r ) = 4 #(i)ab(t, (J r ) + 0(C) -^ 4 n AB ( asym ) at spatial infinity, 

A 9(l)AB(T , 0" r ) = 4 f]AB(asym) + 4 fyl)AB(l~ , ^) , 

A h(i)AB( T , °" r ) = 0(() — > at spatial infinity, (1.5) 

where ( << 1 is a small dimensionless parameter, the small perturbation 4 /i(i)AB has no 
intrinsic meaning in the bulk and 3 (?(i) rs (r, a r ) = — e 4 g(i) rs (r, cr r ) is the positive-definite 
3- metric on the instantaneous (non-Euclidean) 3-space E r . 

The asymptotic 4-metric allows to define both a flat d'Alambertian □ = d 2 — A and 
a flat Laplacian A = ^ r d% on S r (d A = gfx)- We will also need the flat distribution 
c(a, a) = ^ 5 3 (a, a) = — y^z^y with \a — a\ = a/XL ( aU ~ a ' u ) 2 i wh ere S 3 (a, a) is the 
Dirac delta function on the 3-manifold £ r . 

Therefore we will solve the constraints and the Hamilton-Dirac equations in a fictitious 
Euclidean inertial 3-space identified by the asymptotic Minkowski metric and the solution 
will describe the gravitational field in our well defined Riemannian 3-spaces modulo correc- 
tions of order 0(( 2 ). The instantaneous 3-spaces will deviated from flat Euclidean 3-spaces 
by curvature effects of order 0((), in accord with the equivalence principle. 

We assume that the dimensionless configurational tidal variables Ra in the York canonical 
basis satisfy the following conditions 
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\R- a (r,a) = R (1) - a (T,a)\ = 0(() « 1, 

\d u R- a (r,a)\ ~ io(C), |9 u ^i? s (r,,T)| ~ -^O(C), 

|<9 r it> s | = ^O(C), |^it>,| = ^O(C), |9 r 9 uJ R,| = 0(C), 

=}► g a (r,a) = ^-^(^ = 1 + rW(r,<r) + 0(C 2 ), 

r« = £ Taa R„ £ ^ = > ^ = E 7*. ^ ( L6 ) 

a a a 

where L is a 6i<? enough characteristic length interpretable as the reduced wavelength X/2n of 
the resulting gravitational waves (GW). Therefore the tidal variables R a are slowly varying 
over the length L and times L/c. This also implies that the Riemann tensor 4 Rabcd, 
the Ricci tensor 4 Rab and the scalar 4-curvature 4 R behave as 0(C). Also the intrinsic 
3-curvature scalar of the instantaneous 3-spaces £ r is of order -^O(C). To simplify the 
notation we use R a for R(i) a in the rest of the paper. As shown in paper II, this condition 
defines a weak field approximation. 

Eq.(1.5) can be implemented if we add the following assumptions 

= 6 = y/det3g rs = l + 6(l> (1) + 0(C 2 ), 
N = 1 + n = 1 + ri(i) + 0(C 2 ), ra (o) = n (1)(a) + 0{( 2 ), 

4 5 r (i)rr = e + 4 h {1)TT = e (1 + 2 = e + 0(C), 

9(l)rr = A hd) rr = — efj(i)( r ) = 0(C), 

= -^r S + 4 /i(i)r S = -e [1 + 2 (r« + 2 (1) )] 5 rs = -e5 rs + 0(C), (1.7) 

5 ra (1 + if } + 2 (1) ) + 0(C 2 ), respectively. 

As shown in paper II, with these assumptions we have 1 

871 ° n s (r, a) = ^ n (1)B (r, a) = \ 0(C)= [d T R- a -"£ Taa d a n (1)(a) ] (r, a) + ± 0(C 2 ), 

a 

<7(o)(a) = ff (l)(a)(a) = E 7 « a n « a + J 
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while the triads and cotriads become 3 e r mr a ) — K (1 — — 2(f)(1)) + 0(C 2 ) and 3 e(x)( a ) r = 



c 3 



a 



1 Let us remark that everywhere ^-(i)a appears in the combination ^ II(i)a — x ^(0, which behaves like 
d T Ra, i.e. it varies slowly over L. 
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where C( a )( a ) are the diagonal elements of the shear C( a )(&) of the congruence of Eulerian 
observers, whose 4-velocity is the unit normal to the 3-spaces S r as 3-sub-manifolds of 
space-time (see Subsection IID of paper I). For the non-diagonal elements of the shear, for 
the momenta n\ and for the extrinsic curvature the assumptions of paper II are 



0(a){b)\c#b = <7(l)(o)(6)|o#> = L O (0, 



5 X 



12ttG 



7T7 
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(i)« = ^ °(0 



= (1 - <5 rs ) <T(i) (r)(8) + 5 rs 3 iT(i) - d T T« + £ ($ra ~ ^) 5 a n ( i )(a) ] + j 0(C 2 ). 

a 



;i.9) 



Let us now consider our matter, i.e. positive-energy scalar particles described by the 3- 
coordinates r/[(r), % = 1, .., N, such that their world-lines are xf(r) = ^ M (r, ffi(r)). Ki r (r) are 

the canonically conjugate 3-momenta. We have ?7[(t) = 0(1) and t)[(t) = ^h^L d = -l^I = 

0(1) since r = ct (in the non-relativistic limit we have ff i = Vi/c = 0(1) — > c ^oo 0). 

As shown in paper II, to get a consistent approximation we must introduce a ultraviolet 
cutoff' M on the masses and momenta of the particles so that the mass density and the mass 
current density (see the next Section for the energy-momentum of the particles) satisfy the 
following requirements 



M(r,a) = M^ V \r,a)+n {2) (r,a), 



mi = MO{(), 



d 3 aM^ V \r,a) = McO(0, 



d 3 aTZ (2 ){T,a) = McO(( 



M r (r,a) = M {1)r (r,a), J d 3 a M^ V r \r,a) = McO((). (1. 

Here M is a finite mass defining the ultraviolet cutoff: M c 2 gives an estimate of the weak 
ADM energy of the 3-universe contained in the instantaneous 3-spaces E T . The associated 
length scale is the gravitational radius Rm = 2M 4 ~ 1CT 29 M. 



6 



The description of particles in our approximation will be reliable only if their masses 
and momenta are less of McO(() and at distances r from the particles satisfying r = 
| <7 — fji(r)\ » Rm (that is at each instant we must enclose each particle in a sphere of 
radius Rm and our approximation is not reliable inside these spheres). 

Therefore for the particles the validity of the weak field approximation requires 

m(r) = o(i), ^ = 0(1), ^ = 0(0, ^<o(C). (i.ii) 

rriiC Mc M 

Our results are equivalent to a re-summation of the post-Newtonian expansions valid for 
small rest masses still having relativistic velocities (^^ = 0(1), ^ = 0(1)). 

Since, as said in Subsection HE of paper I, we have that the matter energy-momentum 
tensor satisfies Va T ab (t, a) = due to the Bianchi identities and since A gAB = 4 VAB(asym) + 
0((), we must have 8a T^(t, a) = + 8a T^f 2 ) ■ At the lowest order this implies 

8 T M { KP + 8 r Mf^ ] = + 8 A Ufa, 

8 T M { $? + 8 S T(7) = + 8 A Kfa, (1.12) 

as in inertial frames in Minkowski space-time. The equation 8aT£?(t,g) = + dATZf 2 f 
implies 8a (t^ b (t, a) a c — T^(r, a) cr B ^j = + 8a C (angular momentum conserva- 
tion) . 



In conclusion, since the weak field linearized solution can be trusted only at distances 
d » Rm from the particles, the GW's described by our linearization must have a wavelength 
satisfying A ~ L > d » Rm (with the weak field approximation we have A << 4 TZ without 
the slow motion assumption). 

If all the particles are contained in a compact set of radius l c (the source), the frequency 
v — j- of the emitted GW's will be of the order of the typical frequency u s of the motion 
inside the source, where the typical velocities are of the order v ~ u s l c . As a consequence 
we get v = j m u s « v/l c or A ^ ^ l c » R M , so that we get f ^ l f « and l c » R M 
iff =0(1). 

If the velocities of the particles become non-relativistic, i.e. in the slow motion regime 
with v « c (for binary systems with total mass m and held together by weak gravitational 

forces we have also - c ps \J~^ « 1)> we have A >> l c and we can have l c ^ R M . 

As shown in paper II, this HPM linearization allows to get a consistent description of 
GW's in non-harmonic 3-orthogonal gauges reproducing their known properties in harmonic 
gauges. 
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C. Outline of the Paper 



In this paper we look in detail at the properties of the PM space-times identified by our 
HPM solution and we will study the equations of motion of the particles. It will be shown 
how all the relevant gravitational quantities depend upon the York time, which is the general 
relativistic remnant of the special relativistic gauge freedom in clock synchronization. It will 
turn out that they (with the only exception of the ADM Lorentz generators) depend upon 
the gradients of the spatially non-local function 3 /C(i) (r, a) = ^ 3 i^(i) (r, a) (the non-local 
York time) of the lowest order component 3 if (i) (t, <j) of the York time. This study will be 
done in our family of 3-orthogonal gauges, where the Riemannian instantaneous 3-spaces 
S T have a diagonal 3-metric but still depend on the arbitrary numerical function F(r, a) 
determining the inertial gauge variable 3 K^(r, a). 

We will determine the explicit dependence of the proper time of time-like observers, of 
the time-like and null geodesies, of the redshift of light and of the luminosity distance upon 
the York time in these PM space-times. 

Then we will study the consequences of the HPM linearization on the equations of motion 
for the particles and we will make their Post-Newtonian (PN) expansion at all |P7V orders. 
In particular we will show that at the astrophysical level there is a 0.5PN contribution 
to dark matter coming from the relativistic inertial effect connected to the choice of the 
function 3 /C(i). In the non-Euclidean 3-space there is an effective inertial mass different 
from the gravitational mass, even if the equality on inertial and gravitational masses holds 
in the 4-dimensional space-time in accord with the equivalence principle. 

In Section II we rewrite the solutions of the constraints and of the equations for the 
lapse and shift functions (all determined by elliptic equations inside the 3-spaces E T ) in 
3-orthogonal gauges, given in paper II, restricting them to the case in which the matter 
consists only of point particles. An equal time development of the retarded solution for the 
GW's is also given. Then, after the expression of the ADM Poincare' generators we give 
the Christoffel symbols and the Riemann and Weyl tensors of the PM space-times. Finally 
the expression of the 4-spin and 3-spin connections and the expression of the Ashtekar's 
variables in the York canonical basis are given. 

In Section III we show properties of the PM space-times and of their Riemannian 3- 
spaces in 3-orthogonal gauges: the proper time of time-like observers; the 3-distance, the 
3-curvature and the extrinsic curvature of 3-spaces; the PM time-like 4-geodesics. Also the 
comparison of the HPM solution for the 4-metric in 3-orthogonal gauges with the astro- 
nomical conventions for the 4-metric of the Solar System in a suitable harmonic gauge is 
given. 

Section IV is devoted to the PM null geodesies, the red-shift, the geodesic deviation 
equation and the luminosity distance of PM space-times. 

In Section V we give the PM equations of motion for the particles and discuss their 
qualitative properties and differences from the case of charged particles plus the electro- 
magnetic field in Minkowski space-time. After a discussion of the problem of the center 
of mass and of the relative variables in the PM space-times (using the two-body problem 
as an example), we study the PN expansion of the equations of motion and we compare a 
1PN binary system in 3-orthogonal gauges with the standard one in harmonic gauges. In 
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Appendix A there is the 1PN expression of the ADM Poincare' generators with terms till 
order 0(( 2 ) included. 

In Section VI we show that the non-local York time allows to describe the dark matter 
present in the masses of galaxies and their rotation curves as a relativistic inertial effect 
absent in Newton gravity. 

In the Conclusions we comment on the gauge problem in GR, discussed in the Conclusions 
of paper II, and on how our explanation of dark matter may help in finding the PM extension 
of the Celestial reference frame (ICRS) [5] outside the Solar System. Then we show which 
lines of development are opened by our formulation, especially in cosmology and in particular 
for the possibility to reformulate the back-reaction approach for the elimination of dark 
energy, based on averages of scalar quantities in 3- volumes of 3-space, in the York canonical 
basis. 
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II. THE PM SOLUTION FOR THE GRAVITATIONAL FIELD 



In this Section we review the results of paper II when the matter consists only of point 
particles. At this order the HPM linearized solution in the family of 3-orthogonal gauges 
depends on the York time 3 K only through the following function 3 /C = ^ 3 K. Only the 

linearized ADM Lorentz generators have also a dependence on 3 K. We give also the equal 
time development of the retarded solution for GW's. The metric, Christoffel symbols, spin 
connection, Riemann and Weyl tensors of PM space-times are given. Finally Ashtekar's 
variables are expressed in the York canonical basis and their PM limit is found. 

A. The Energy-Momentum of the Particles 

From Eqs.(3.9), (3.10) and (3.12) of paper II we get the following expression for the 
energy-momentum of the particles (in the following equations the notation -p- 0(( 2 ) means 

Ef =1 5\aMr))MO{e)) 

N 



M^ V \r,a) = 7ft(r, S) = E 5*(a, #(r)) Vi yfij + ««(t) 



i=i 

N 



£V(ff,i*(T))ift g| +§Q(C 2 ), 



N N 



M (1) 



i=i i=i Jl- n i (r) 



w — V a / m.. c -4- k.. t / 



\/ m f c2 + ^( r ) 



W,^, u ^0 (1) C(r) + E a tt 1 ' fa a (r)) 2 (r) V 
= 5 (a,^(r))^m iC (^ J(t,o-), 

<=i ' V 1 — »7i (t) 



= ^(^)^5 3 (M l (r)) J)t « ir (r) = E ^(^ ^ -^^L + MO(( 2 ), 



<=i i=i y 1 - ^ (r) 

AT AT r/ . v . s( v 

i=i i=i Jl _ 7) ( r ) 
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N 



K ir(T) K is\J) 

\Jm\c 2 + k 2 (t) 



A? 



\rs 



Ki r [T ) Ki s [T ) 

yjm 2 c 2 + k 2 (t) 



1=1 



K 2 ( T ) 

A^W (T ' a) = -^^ 4vr|a-^(r)| ' A E T « (r ' a) = " E 4. |<? - ffr)\ ' 

AT 



K,v ( T } 



I E d c M^ V c \r,a) = "E * E Mr) / (^Iff-^.^^is 

3 K-#))K-#)) 
l«?i-^(^)l 2 



(2.1) 



where we used /tj = + McO{() and A/m 2 c 2 + «? = + McO(C). We have also 

given the second order of the mass function. The quantities g' TT , g r ' TS and g' rs are the mass 
monopole, the momentum dipole and the stress tensor monopole respectively (see Appendix 
B of paper II). 



B. The Solution of the Super-Hamiltonian and Super-Momentum Constraints and 
the Lapse and Shift Functions for the Family of 3-Orthogonal Gauges 

From Eqs.(4.6), (4.7), (4.16) and (4.17) of paper II we get the following expressions for the 
solutions: a) 0(i)(r, <r) of the super-Hamiltonian constraint; b) N(r,a) = 1 + n^{r,a) and 
n(i)( a )(r, a) of the equations for the lapse and shift functions in the family of 3-orthogonal 
gauges; c) o"(i)( a )(&) \ a ^b(j, er) (the off-diagonal terms of the shear of the congruence of Eulerian 

observers) of the super-momentum constraints (see Eq.(1.9) for rcf^): 
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4>(r,a) = 1 + 60 (1) (t )( t) + O(C 2 ), 



0(i) (t, a) 



2 " ? >s v, +ie!^1(^)= 



c 3 A 



G 



A/mfc 2 + /5f (r) 



a (Ti — 



<7- W"J 



16tt 



167T J \<J — G\ 

d a i zr-j 



(2.2) 



7l ( i) (T, <?) 



4ttG 1 



d T 3 JC 



(i) 



c 3 4^ * |<T-r7i(r)| 



( 1 + s^)-* ,K <" (T ^ 



(2.3) 



= [^ 3 %)+ 1 ^i(4<r ) -|E^< ) r ) )+ 

c 

+ ^|^(4r«-$:|r«)](r,a)= 

c 

l^-^r))^)-^-^))- 



\v-Vi(T)\ 2 



<i 3 (Ti 



47T |(T — (7 i 



<9la 



E c 5 2 2 c r«(r,a 2 ) 



87T C?i — C?2 
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= 9 a 3 /C (1) (r,a)-^— * 



rrii 



l(a a -r 1 t(r))ff i (T)-(a-ff i (r)) 



d 3 



471 la — ai\ 



9lad T 



2r< 1 >(r,* 1 )- /dV 2 



, E c 9 2 2 c r«(r,a 2 ) 



87T C?i — (72 



(2.4) 



0" (l)(a)(6) |a#>(T , <?) = ^ ^WOO + ^ b n (!)( a )) l«^( T ' ^) = 



(1) 



+ s^(r« + r«-ij:|r«)](r, 3 )= 



1 <9 a (9fc ,({7V)" 

2 A ^ A « c . 

c 



+ 



l(^-^(r))^(r)-(<T-r/,(r)) 



|0 - Vi{T) 



+ 



+ 



d 3 



(J\ 



4-7T \a — a-A 



d 



id 



(i) 



2d T T d 



8tt l&i — &2\ 



+ d a d b 3 )C il) ( y T, a). 

(2.5) 



The act ion- at- a- distance part of the solution is explicitly shown. Only the PM volume 
element 0(i) = 1 + 6(f)(1) is independent from the York time. Eq.(2.4) describes gravito- 
magnetism in these PM space-times 2 : it has an inertial gauge part <9 a 3 /C(i). 

C. The HPM Gravitational Waves 

By using Eqs.(7.1) and (7.2) of paper II, the retarded solution for the tidal variables 
with the matter restricted to point particles is (see Eq. (3.12) of paper II for T rs ; the TT 



2 The gravito-magnetic potential Aq has the components Aqi t \ ~ c 2 ni\\i r \. The gravito-magnetic field 
Bo(r) — c ^G(r) = (d x Ac) r is proportional to the second term in the Christoffel symbol r" 1 ^ )rr given 
in Eq.(2.15). Instead the gravito-electric potential is $ G = -^n {1) = -^f- ± {Mf x p + J2 a T (Tp + 
^- d T 3 /C(!)). Both A G ( r ) and $ G depend on the non-local York time. 
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projector V rS uv is defined in Eqs.(6.7) of paper II; the spatial operator M a ^ is defined in 
Eqs. (6.24) and (6.25) of paper II) 3 

a 

2G . ~ 



C 2 



b 


i 

v-*i\)fli 


uz 

[r- 


a - (Ti|) 5 3 (ai - ffi(T - | 








a — u\ 





J 1 Jl-tfir-W-a^) \3-Bi\ 

(2.6) 



By making a equal time development of the retarded kernel like in Ref . [7] for the extrac- 
tion of the Darwin potential from the Lienard-Wiechert solution (see Eqs. (5.1)-(5.21) 
of Ref. [7] with £ s P r ± s (a)f)!(T) ^ J2uv Vbbuv{o) ^ffffi ) and by using the fact that 

fji(r) = 0(C) (see also Section V) we get the following expression of the HPM GW's from 
point masses 



r<"(r,if) 



C b i uv V 1 _ ^ l r J 



9a 



+ 0(C 2 ).(2.7) 



As shown in paper II the multipolar expansion of the TT 3-metric and of the tidal 
variables reproduces the quadrupolar emission formula (q uv \ TT is the mass quadrupole; the 
TT projector A abcd is defined in Eqs. (7.17) of paper II) 



3 One could study the radiative fields r„ (r, a) at null infinity (|<?| — > oo with the retarded time r — \a\ 
fixed) to see whether terms in In \a\ appear like in the standard approach to GW's in harmonic gauges 
(see Section 5.3.4 of Ref. [6]), but this will done elsewhere. 
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1 ^0> S ( T > °) = ~ e E A rsuv(n) UtH — + (higher multipoles) + 0(l/r 2 ), 



02 ? u W |tt 


((r- 






— * 





^(r.o) = -i£^^V) ^ r """|V V " — + (higher mulUpoles) + 0(l/r 2 ), 



aft 



t—iuv bbuv U T H 




*l) 




a 





^ / 

= E * ^ r - \*\) ^ T - i*d \K c2 + - i*i) = 



i=l 
IV 



\- mjCrif(T - |<?|)^(r - 

= E ^ 1 , 2 • (2- 

i=i yi-rj^T- \a\) 

Eq.(4.18) of paper II gives the following expression for the tidal momenta of Eqs.(1.8) 
(namely for diagonal elements <7(i)( a )(a) of the shear, see Eq.(1.8)) 



8tt G 



n B (r, <?) = 9 r R- a (r, a)-J2 Taa \d T ^ (4 r« - 1 £ <9 2 lf>) + 



c 



c 



c 3 A ° ^ a A 



= (E - E i - ± E »*m%?) + 

ft a c 

+ d 2 a 3 IC {1) ])(r,a), 

M- al = 5- ai ~J2 Taa f (2 Tba ~ 3 E T» ^) , (2-9) 

a ft 

While the tidal variables R & do not depend on the York time, the tidal momenta n s 
depend upon it. 

D. The Weak ADM Poincare' Generators 

The HPM linearization of the ADM Poincare' generators is given in Eqs. (4-21)-(4.24) of 



paper II. The final expression of the generators is (\Jmfc 2 + k 2 (t) = c/y 1 — %(t); the 
spatial operator M sb - is given in Eq.(2.6); L is the GW wavelength of Section III of paper II) 
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M (1) c + - E A dm(2) + McO(( 3 ) = 



- £ 



K7 



(r) %^§HP+£ C ^r«(r^(r))l +E C lf>(r, *(r)) 



v/mfc 2 + /tf(r) 



G 



E M 



A/m 2 c 2 + /t 2 (r) ^m 2 c 2 + k 2 (t) 



1^0") -^(l") 



4^(r) ■ Kj(r) 



— > K ir {T) K js (T) / d (J — — =T7 x.o + 



E m H-r) ■ 9 3 /C (1) (t,7t;(t)) + McO(C 3 ), 



(2.10) 



4m = f {l) +p\ 2) + McO{e) = 

3 /* 

= J] Vi «ir(r) - ^ J d 3 aJ2 ( d r R* ^ d T i2 6 ) (r, <?) + 

+ E * E ^ ( E ^ r ^ - 2 r i 1} ) *oo) - 

j a c 

i 
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jr? )+j rS )+M cO(( 3 ) = 

Yl (pK T ) K is(r) - (r) «ir(r)) - 

i 

2 E*E«»w ('. r w^-'.'M8?) 

in li li 

|(r<%«W)-iE|r«) + 

c 

2 5> hM |(r?>- 1 £§!*>)- 

i c 

- ^MKrW-i^fr?))]^,^))- 

c 

Jd 3 a[j2 W r d s - * s d r ) R- a M- al d T R~ b + 2 3 K (1) d r d s (lf> - r«) 

r, <?) + 



+ 



8ti G 



+ 



+ 2 (a, r?) + s T rf ) -^^ if>) ^ (rf ) - rp>) 



+ £ % y/m?<? + i%(T) (<(r) A - tf(r) ^) 3 /C (1) (r, #(r)) + McLO(C 3 ), 



9 



(2.12) 
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-J ADM = J(T) + j(2) + MC L0(( 3 ) = 



G v ^( r ) V m j c2 + ^( r ) 77[(r) 

_ _ (ri 1} + \ E c f r«)(r,^(r)) 

+ ?**< T) ? 

-8^E^r<V,.-)^(r«-I E fr«)( Tl(f )- 

a c 
a c d 

+ i E E f «. (r« - i E f rf '(r, *>) + 



c 

+ 



lfaG E(^(^rf>(r.<?)) 2 + 



a,b 

.3 rflftfl / ... 1 



+ I^E[^(rl" + r«-iE|r-)(^-)] 

1 ^ d a d b d T f (1) (1) 1 ^ <9 c 2 lf ] x ^ Kia ( r ) (<x a - 7/f (t)) 



a^b c 
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c 3 



0,6 i 1 

1 V- V- ^ b (r) (a a - 77? (r)) <9 a <9 6 / 3 G ^ g^r) ■ (g-^(r)) \ 

V- ^(r) ■ Kj(r) (g - 7fj(r)) ■ (g - g(r)) + ^(r) ■ (g - ^(r)) ^(r) ■ (g - ^(r)) 1 
+ 2vrc 3 ^ ^ I ^ _ ^.( r ) |3| _ ^.( r ) |3 J 



+ / «■» [| E *^$f - i^ 1 ' (E f^OM - 

" / ^* [ijfeK 1 ^)' " E (rS 1, (^)' " l(E f ri»M) 2 ] - 

s c 

G v/m^c 2 + ^(T)Jm 2 jC 2 + /? 2 (r) 

(2.13) 

At the lowest order they reduce to the special relativistic internal Poincare' gen- 
erators in the rest-frame instant form of Ref.[4]. They are = M^c = 

i(T) = Ei ViVi{ T ) \A«iC 2 + k 2 (t) rs 0. The conditions j'^j « and « are the rest- 
frame conditions eliminating the 3-center of mass and its conjugate 3-momentum inside the 
3-spaces of the rest frame. As shown in Ref.[4] in special relativity (and also in PM canonical 
gravity) there is a decoupled external (canonical but not covariant) 4-center of mass to be 
used as collective variable. 

Eqs.(2.12) and (2.13) show that the Lorentz generators depend both on the local and 
non-local York times at the second order. 

For the effective Hamiltonian for 3-orthogonal gauges, not used in this paper, see Eq.(4.26) 
of paper II. 
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E. The 4-Metric, the Triads and Cotriads, the S T - Adapted Tetrads and Cotetrads 



Eqs. (1.5) and (1.7) imply the following expression for triads, cotriads, tetrads, cotetrads 
and the 4-metric [lf o) = (1; 0); el( ) A = (1; 0)) 



e (l)(«) 


= ^(i-rW-20 


(i))' 


3 e ( l)(a)r 


= 5ar(l + r« 


o A 








o A 


4 F 


= ^) + l ti) = (i - 






4 F 

^(l)(o) " 


o(o) 








o (a) 


4 F 


= e (^(o)A + ^(i)a) = 


(1 + 


n(i))(l;0) 


4 P 

> ^(1)^ — 



e 9(i)tt 



9(1)7 



-£ 9{l)r 



n 



(i)( 



•3p ^ 

a) i e (a)r I 5 



l + 2n (i) = 1 



8n G 1 



~ 3 A 



•^S y) +E T ( a i))- 2 ^ 3 %)' 



<9 r 



-"(l)(r) = 

- — 9 ^4rw - V ^r^A 

2 A V A c /' 

c 

3 ^ ( i)r S = ^ [i + 2(r( 1 ) + 2 (1) )] 



1 - 



c s A W r ^ A 



(2.14) 



o A o_A 

The tetrads i E^ Vj ^ are adapted to the 3-spaces: 4 -E(i)( D ) = lf \ + i s the normal to S T . 
While the triads and the 3-metric in E T are independent from the non-local York time, the 
4-metric components A g(\) Tr , A 9(i) T r and the tetrads depend upon it. 



F. The PM 4-Christoffel Symbols and the PM 4-Riemann and 4-Weyl Tensors 

By using the PM linearized 4-metric given in Eq.(2.14) we can evaluate the Christoffel 
symbols and the Riemann and Weyl tensors of these PM space-times and study the properties 
of the Riemannian instantaneous 3-spaces. While the terms containing M^ly \ T^, 

correspond to act ion- at- a- distance contributions, the terms containing = 'jar R& de- 
note retarded GW contributions. The non-fixed gauge part is given by the terms depending 
upon 3 £ ( i) = ^ 3 ^(i). 
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1. The PM Christoffel Symbols 



The Christoffel symbols and their linearization have the following expressions in our 
gauges 4 ( 3 r" s is the Christoffel symbol built with the 3-metric 3 g rs of the 3-space) 

^ A bc = \ V* (d B 4 9ce + d c 4 g B E - d E 4 g BC ) = 

= Vdbc + 0{C) = \ V s (d B 4 g { i)cE + d c 4 g { i )B E - d E 4 g {1)B c) + 0(C 2 ), 



4pT 

TT 



(pr n + n (a) 3 e\ a) d r n - n (a) n {b) 3 e s (a) 3 K rs 3 e[ fe) ) = 

4 r[ 1)TT + o(c 2 ) = d T n (1) = ^ ±d r (m^P + - dl 3 /C (1) + 0(C 2 ), 

a 

A Kr = -^{d r n- 3 K rs 3 el a) n [a ^ = 

A n 1)rr + o(c 2 ) = d r n (1) + ^ | (-mS^ + E T «) - * 3/c « + °(c 2 )^ 



4 r T = — 3 k = 

l + n rs 



(1) 

AttG 1 / . jttv\ „ . .rr/v^ ct <9. 



4 r[ 1)rs + 0(C 2 ) = (9 r n (1 ) W + 9 S n (1)(r) ) + 5 rs d T (r« + 2 + 0(C 2 ) 



: 3 A 

+ 5 T r?) - ^ d T (r?) + rf )) + 1 (5 rs + ^) s T £ g r« - 



c 

d r d s 3 !C {1) +0(( 2 ), 



4 For the evaluation of 4 r" r we need the kinematical Hamilton equations for the cotriads given in 
Ref.[8]. If we are in a Schwinger time-gauge with ar a \(r,a) 7^ 0, we have to add to 4 r" T the term 
~\ J2abc 3e ('a) e (a)(b)(c) A(6) n (c), where fr( a )(T, a)=d T (*(«,) (r, ct) is the Dirac multiplier in front of the rota- 
tion primary constraint. 
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d T n 



Kt = 3e h dr n («) ~ n(a) + ^ n(h) n(c) 3 ^ a) 3e " fe) ^ s 3e(c)r ~ dr 3e{c)s " > + 



bcrs 



n>(a) ri(b) 



+ E %) H) + (1 + n) %) (*(•)(*> (1 + n) 2 



<9 r n 



bcrs 
1 



rb 



= 4 rf 1)rr + 0(C 2 ) = d T n (1)(u) + <9 U n (1) + 0(C 2 ) = 



(i) 
47rG 



I - x E + x ("So" + E iS) 



<9« 



3» 



A 



A 



+ J^W-5]f r«) + o( C 2 ), 



4p« 3^m 

1 rr — e (a) 



■n (a) - (1 + n) (0(a)(6) - 



1 + n 



(a.) u (b) \3t^ 3-s i 3 s 3 C o _ 



(1 + n) 2 

+ 3 c V (a) 3 ^(c)r 3 e( b) (d v 3 e (c)s - <9 S 3 e( c)t) ) j n (o ) = 

= 4 r^ 1)rr + 0(C) = 5 ur d T (r« + 2 (1) ) + l - (d r n (1)(tt) - 9 tt n (1)(r) ) + 0(( 2 ) 

^r^l(i)( u ) C '«- /Vt (l)(r) 2 UT T (!) / 



c 3 A 



+ 5 ur d T (r?) + (r?) - rW) + o( C 2 ), 



4 "pit _ 3pw , n 0) 3 s u 3tz _ 
1 rs ~ L rs + (a) ~~ 

= 4 r^ + o(c 2 ) = 3 r^ 1)rs + o(c 2 ) 



<W ^ (r« + 2 + 5 US 9 r (r« + 2 - 5 rs a (r« + 2 </> (1) ) + o(c 2 ) 

_ Air G 5 ur d s + 5 US d r - S rs d u .auv) 



A 



+ (s ur d s + 5 US d r ) r« - 5 rs d u r« + 6 ^9s + Su S d r 5 rs d u ^ ^ ^ + 0(c ^ 



(2.15) 



U. T/ie PM Riemann and Ricci Tensors 

The 4-Riemann tensor and its linearization have the following expressions 
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J^ABCD — 9AE ti BCD — 

= ~\(pA d c 4 9bd + d B d D 4 g AC - d A d D 4 g BC - d B d c 4 fiUr>) + 

-r 9EF I J- AD 1 BC ~ 1 AC 1 BDI ~ 

= A R(1)ABCD + 0(C) = A VAE 4 Rfi)BCD + 0(C) = 

= ~ (d A d c 4 g(i)BD + d B d D A g( 1)AC - d A d D 4 g(i )BC - d B d c 4 #(i)ad) + 0(C), 



4p —t 3 Ft 



— e 



5 rv d s d u (if > + 2 - 5 ru d s d v (if > + 2 + 



+ 5 SU a, d v (if) + 2 (1) ) - <y w a, a, (if> + 2 (1) ) 



4tt G (5 rv d u - 5 ru d v )d s - (5 s v d u - 5 su d v )d r j^uv) + 



A 

+ (<U a* - 5 ru d v ) d s if) - (<j w d tt - <5 SU a,) a, if) + 



+ 



2A 



4p _ , 



(5 rv d u - 5 ru d v ) d T (if) + 2 0(i)) + 2 d r ( 9 « 

^(l)(w) 

^ ^ ((5™ ft - 5™ ft) d r - 2 ft (ft jug? - ft M%£j) + 

+ ft ((5 rv ft - s ru ft) (if> + 1 £ f if >) + ^ (if) - if))); , 



= e 



4n 

- rl '(l)T7'TS 



= e 



+ ft 



-| (2 ft ft n (1) - 2 <5 rs ^ (if > + 2 + ft (ft n (1)W + ft n (1)( , 
^ i (ft ft + ^ T^) + dl MfP + 



2(ft*C +4^) + 



+ a 2 (*„ (if > + i £ § if >) - ^ (if> + if) - i £ § if>)) 



= e 

+ a 
1 



4ttG 1 



2 3 A 



2 (a, M ™ + ft - ^ £ ft ) - 



2 °V "'(l)rs- 



(2.16) 
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The final expression of A R^ Trrs has been obtained by using Eq.(6.12) of paper II and has 
been used in Eq.(7.39) of paper II. Let us remark that the Riemann tensor does not depend 
upon the York time 3 K. 

For the 4-Ricci tensor and the 4-curvature scalar we have (□ = d 2 — A) 

4 R(1)AB = 4 f] EF 4 R(1)EAFB = £ (^R(1)tAtB ~ 4 R(l)rArB j , 

r 

4 R(1) = 4 V AB 4 R(1)AB = e (^R(1)tt - E 4 -^(l)rr), 



4 J R ( i )rr = -6 d 2 T 0(i) + A n ( i) + E d T d r 



V n (l)(r) 



c 3 



** G ((1 + 3 f ) M ™ + E T$ + 3 | E *e *C) , 



i? ( i )rr = <9 T <9 r (r^ - 4 0(i)) + i E d s (^(i)(.) - 9 s n(i)( r )) = 

s 



3 V A « 



4 



i2 ( i)„ = <9r d s (-ri(i) + lf> + If) - 2 0(d) + 5 rs (d 2 T - A) (r« + 2 

- ^ T (9 r n ( i Ks )+9 s fi(i)( r) ) = 

- A U % TT \ AnG ( 5 -M {uv) u rr 

a 

- 2 1 (s rM ^j + a xsr») + e a c Mgr») , 



'Ad, = 2 ( - yj a? r<" + a n( i) + yj a x a r n (1)w + 8 a <a (1) - 12 a; <a (1) ) 

r r 

.(rrv) 



c 3 



(2.17) 



By using Eqs.(1.12), it can be checked that Einstein equations 4 Rab~ \ 4 9ab 4 R=^f- Tab 
are verified, namely we have 4 R(i )A b ~ \ 4 Vab 4 R(i) = ^T {1)ab + 0(( 2 ). 

3. The PM Weyl Tensor 

For the Weyl tensor and its electric and magnetic components with respect to the Eulerian 
observers, whose unit 4-velocity l A is the normal to the 3-space S r with the zeroth order 
expression = (7^1;/^ = 0) [see Eqs.(2.14)], we have 
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A Cabcd — A Rabcd — 2 ( A 9ac 4 Rbd + A 9bd A Rac — A Qad A Rbc — A 9bc A Rad) + 
+ g ( A 9ac A 9bd - A 9ad A 9bc) A R = A C(i)abcd + 0(( 2 ), 

4/~i 4/~i 4/~i _ 4/~i 

t-'ABCD — t-'CDAB — ~ ^BACD — — ^ABDCi 

a Cabcd + a Cadbc + a Cacdb = 0, 



i C(l)rrrs — A R(l)rrrs ~ ^ R(l)rs ~ $rs A R(1)tt j — g ^rs 4 -R(l) — 



4:1 TT 



"4 (° " A ) 



4ttG 



_ 4p _|_ _ | A 4 B -A 4 R 

^(l)rriti; - n '(l)Truu ' I w nj - ri '(l)TU "rti -^(l)-™ 



= ct 



(<U 5 M - s ru d v ) (r?) + ± £ § + ^ - r ^)] + 



+ 



C(l)rsM-i) R(l)rsuv ^ {^ru R(l)sv &sv R(l)ru 



8 rv 4 R(l) su — 8 SU A R(l)rvj + g ^ 



1 nrA 4 h TT -i- a 4 h TT — a 4 /) TT -A 4 /) TT 



1 d 2 

- (S rv d u d s + 5 SU d v d r - 5 ru d v d s - 5 SV d u d r ) (F^ + - £ -J T^) + 



+ 



47rG r,2 ,„ „ „ <9 n <9 s , r d v d r ^ d v d s ^ d u d r r uv) 



A 



A 



A 



A 



, , , ,-, , d u d s d v d r d v d s d u d r ^ d T d c 

+ {o ru S v - O rv O su - 2 d rv — 2 b su — h 2 d ru —— + 2 ——J —J— M 



A 



A 



+ | (S ru (d s M^ + %M<™) + 5 SV (drM™ + d u M%; 



(UVh 



l (l)r 



25 



^(1)5 ~~ *7 ^{Y)CEBF t( ) t( ) — ?7 tv(l)C*rBr — — £ "Bs ^(l)-rrrs, 

rs 

H(i)AB = 7s e AMCD ifo) i f] DE 4r ] CF 4 C(i)EFBG lf ) = 9 ^( 



7 (l)™rs 



X) ^ ^ (d T [e rsu d u (lf> + ^ X) a r ^ )} 



<9 r <9 W <9^ ^ 



A 



4ttG 1 



~ 3 A 



(2.18) 



Their Newtonian limit, in particular the vanishing of H^ab, is consistent with Ref.[9]. 
G. The 4-Spin and 3-Spin Connections 

The 4-spin connection 4 wa^| associated with the general tetrads 4 E^ is connected with 

the E T -adapted one 4 Wyi|^ by means of the Lorentz boosts with parameters </?(„) [8]. The 

expressions of these 4-spin connections are {jl( a ) ( r ) <?) = 9 T a( a ) ( r ) <?) are the Dirac multipliers 
in front of the primary first class rotation constraints; they vanish in the gauges a( )(r, a) « 
0; 3 K rs is given in Eqs.(2.3) of paper II) 



4 , , O) 



4,°, (a) 

(/?) 



B \°A ^(f}) + 1 AC ) ~ 

L((p( a )) 4 u A L _1 (v?( a) ) + d A L{<P(a)) L~ l ((f {a) ) 
o (a) / o S o C" s 

[pA 4 E {f3) + T% 4e (i3)J > 



(a) 



08) > 



4,°, (o). 



^V"'(a) — e ' w r(o)(a) — 3 e( a) <9 r n - 3 e( a) 3 K rs 3 e( 6) 

4 ° (a) 4 ° 

(6) = -£ W r (a)(b) = — £(o)(6)(c) /%) + 



+ J2 ( 3 4) ^ n M - ^(a) ^ + J2 H) 3 4) 9 » 3e W«) " 

~ \ Yl n w [ 3e U ( 3e w ^ 3e ( a ) s ~ 3e w ^ 3e w s ) ~ 3e w H) 3e ( c ) s - 9 « 3e ( c )«) 



4° (o) 4° 4° 

W r l ; (a) = £ w r(o)(a) = ~ C W r ( a )( G ) 

4 ° (a) 4 ° 3 

(fe) — -£ w r(a)(6) — w r(a)(6)- 



3 K 



(2.19) 
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The 3-spin connection 5 3 w r ( a )(6) = 4 w r (a)(6) with 3 w r ( a )(b) depending on triads and cotriads 
and also on the angles ar a \ is connected to the 3-spin connection 3 u) r (a)(b) i n Schwinger time 
gauges by local SO(3) rotations R(a( a -\) (R T = R^ 1 ) and has the expression in gauges 
near the 3-orthogonal ones having # l (r, a) = O^Jt, a) = 0(() ^ (we use V ra {6 1 ) = 
5 ra - J2i e rai^(i) + 0(( 2 ), see before Eqs.(2,9) of paper I) 



3 Wr(a)(fe) = e (a)(6)(c) 3 ^r(c) 



R(oi( e )) u r R (a (e )) + R(a( e )) d r R (a (e )) 



(a)(6) 



r(a) 



4e 



e(o)(6)(c) V(6)(c) 



6c 



1 



)(&)(c) 3e (6) 



<9 r 3 e( c ) M - d u 3 e( c ) r + ^ 3 eJ c ^ 3 e( d ) r d v 3 e^ d y t 



bcu 

= ^(l)r(a) + 0(C 2 ) = 
1 



dv 



6c 



w (l)r(a)(6) = e (a)(6)(c) w (l)r(c)- 



(2.20) 



Once we know the PM 3-spin connection 3 W(i) r ( a )(6) in the 3-orthogonal Schwinger time 
gauges, we can go to non-Schwinger gauges near the 3-orthogonal ones (with (X( a \(T, a) ^ 0, 
ip {a) (r,a) ^ 0, P(t,o) = Ol^a) = 0(C) ^ 0) and find (we use Eqs. (1.9), (2.5) and (2.14) 

and n {a) = J2 b %)(6)(«(c)) n (6) , 3 e[ a) = £ 6 #(a)(&)(a(c)) 3 e r (b) ) 



4° (a) 3 3 

UV = UJ r (a){b) = e (a)(6)(c) ^r(c) 



R(a>( e) )d r R (a (e) ) + -R(a (e) ) 3 a; r i? (a (e) ) 

3,-. nT/„ , \ /•>/ 2 



(a)(6) 



^ (a (e) ) + i?(a (e) ) 3 w (1)r ^ (a (e) ) + 0(C 2 ), 

J (a)(6) 



= [-R(Q!(e))9 r 

4 W r (o) ( a ) = ^ -R( a )( & )(a( e )) 3 i^(i) rs + 0(C 2 ), 
6 

4 ^ (a) (6) = " 5^ e W(6)(c) <9 r a (c) + ^ R {a){d) (a {e) ) R {b){c) (a {e) ) (d c n {1){d} - d d n (1)(c) ) + 

c cd 
+ ^2 {^(b)(c)( a (e)) d c R( a )(d) - R(a)(c) 9 C R( b )(d) ) ™(l)(d) + 0(C 2 ), 

(2.21) 



4 w r (°) (a) = ^i2 (a)(b) (a (e) )&n (1) +0(C 2 ) 



5 It is defined by the vanishing of the generalized covariant derivative acting on both types of indices of the 
triad: d r 3 e" a -, + 3 r" s 3 e( a -, + 3 u; r(a )( b ) 3 e" t) = 3 V r 3 e" a j + 3 w r(a )( 6 ) 3 e" b j = 0, so that the 3-Christoffel symbols 

(see also the last of Eqs. (2. 15)) have the expression 3 r" s = \ 3 e" a ^ ^ 3 e(;,) r 3 aj s ( a )( b ) + 3 e^)s 3 ^r(a)(fc)) — 



3 e( a ) r d s 3 e^. + 3 e( a ) s d r ej a 
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FromEqs. (3.5) we have 3 K( 1)rs = <J(i)( r )(s) \r^s+S rs (§ 3 K-d T T { r~ ) +d r n {1) (r)-Y.a 9 a n^ a )j 
with <T(i)t r )(s)\r^s given in Eq. (2.5) after the solution of the constraints. 

For the relation between the 4-field strength and the 4-curvature tensors and between 
their 3-dimensional analogues we have 

4 n (a) 4t-.(7) 4 jp(6) 4 (a) 4pC 4 771(a) 4 tjD 

**AB \/3) = & A L B W '(/3)(j)(8) = R DAB & C h (f3) = 

4 (d) o 4, , (a) , 4, , (a) 4, , (7) 4, , (a) 4, , (7) 

— OA (/3) ~ B WA (J3) + Wa v ; ( 7 ) U B y " (J3) ~ ^B (7) ^A (13), 

4 pi _ 4 rpA 4 771(5) 4n (7) 

ti BCD — ^( 7 ) -&b "CD (5)) 

^rs(a) = 2 6 ( a )( fc )( c ) 3 ^rs(b)(c) = ^ € (a)(b)(c) 3e (b)t ^^c) ' ' ^ 'wrs = 
= 8 r 3 U s ( a ) - 8 S 3 UJ r (a) - e (a)(6)(c) 3 ^r(b) 3 U S (c), 

R T stw = 6(a)(6)(c) 3 e( a ) 5(6) (n) 3 ^ 3 ^tw(c), 3 Rrsuv = 6(a)(6)(c) 3 e( a )r 3 e(f,) s 3 fi ul ,( c ). 



3 nf 



(2.22) 

The first Bianchi identity 3 R t rS u + 3 R t sur + 3 R t urs = implies the cyclic identity 

"rs(a) C( a ) _ U. 

H. The Ashtekar Variables in the York Canonical Basis and their PM Limit 

As shown in Ref. [8, 10], the canonical basis 3 e^ a y, 71"^, formed by the cotriads on S r 
and by their conjugate momenta (see Eq.(2.8) of paper I) can be replaced by the following 
canonical basis of Ashtekar's variables (7 is the Immirzi parameter; we use the conventions 
of Ref. [11]; Q a = e £« 



%) = 3 e 3 e[ a) = 2 / 3 ^i? (a)(fe) (a (e) )V;a(^)e-^^^, 

b 

3 A(~f)(a)r = 3 U r (a) + l 3 e S ( a ) 3 K rs , (2.23) 

with 3 w r ( a ) of Eqs.(2.20) and with 3 K rs of Eq.(2.3) of paper II. This formalism is usually 
defined in the Schwinger time gauges ip( a ) (r, a) ~ of ADM tetrad gravity. 

In Ref. [11] it is shown that we have {A^ a y(r, a), A^^ s (t, a)} = due to the re- 
sults { 3 K (a)r (T,a), 3 uj s{b) (T,a l )} = { 3 X (6)s (T,oV^r(a)(T,^)} = ( 3 A" (a)r = 3 K rs 3 e s (a) ), 
which are a consequence of the fact that 3 e^ <5 3 w r ( a ) is a pure divergence (this implies 
3 w r(a) (r,<r) = [ 3 e {b)r 3 cj (b)(a) ](T,a) = J d 3 ^ Y. b 3 ^(b)(b)(r, a)). We also have 

(a)V ' I 

{ 3 A (7)(a)r (r, a), 3 S s {b) (r, S^} = 7 5 s r 5 {a){b) 5 3 (a - a). 

The SO(3) connection A^ a y is considered as a SU(2) connection with field strength 

3 -^(7)(a)rs = <9r 3 ^( 7 )(a)s — d s 3 A( 7 )( a ) r + e( a )( b )( c ) 3 A( 7 )( 6 ) r 3 A( 7 )( C ) S . Instead the true 
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S0(3) connection, associated with the 0(3) subgroup of the Lorentz group 0(3,1), is 



2 e (a)(b)(c) 



R(a {e) ) d r R J \a {e) ) + R mm) (a {e) ) 3 u r{m){n) RT m)(c) (a {e) )) . 

J (b)lc) / 



(b)(c) 
3cr 



Instead the densitized triad £L\ is considered an analogue of an electric field. 

In the Ashstekar formalism the non-abelianized rotation constraint 3 M( a )(r, a) ~ 
of Ref. [3], whose Abelianized form in the York canonical basis is irW = 
— Eb 3 M(b) A(j,)( )(o!( e )) ~ 0, is replaced by the Gauss law constraint G( a ) = ^ r d r 3 £^ + 
" 1 ,i - C/ - 0. The super-Hamiltonian constraint "H(t, a) takes the 



form 



E 



(7)(o)rs 



;!+7 2 ) Efec e (a)(6)(c) 3 ^ ; 



)(6)(c) "-f^r(b) ~K S ( C ) 



E3cr 3cs I 
mn e (a)(m)(n) £( m ) £( n ) i" 



7 1 (1 + 7 2 ) ^ ar G( a ) <9 r ( 3 e) 2 3 £J\ ps 0, while the super-momentum constraints become 



» 

7" 1 Ea E s 3 *< 7 )<.)r< 3 ^ a) - (1 + 7*) 3 K r(a) G {a) 



0. 



The linearized Ashtekar variables in gauges near the 3-orthogonal ones (with a( )(r, (?) 7^ 
0, na) (r,a) ± 0, ^(r,a) = e\ x) {r,a) = O(C) ^ 0) are 



E U = £ flwmHo) ki-r r (1) + 4 (1) )5 rb -^ er6 ,^J +o(c 2 ), 



(7)(a)r 



9 



)(&)(c) 



R(a( e ))d r R (a (e )) 



6c 



+ 



+ ^ ^(b)(m)(«(e)) -R(c)(n)(«(e)) (L <9 n ~ <W <9m) + 2 0(1)) ~ ^ Cmm 9 r 9\ 



(1) 



+ 



+ 7 5^ R (a){b){a(e), 



a (l)(r)(b)\r 



+b + S rb Q 3 if - <9 r + <9 r n ( i )(r) - d a n ( i) (o )^ ) . 

a 

(2.24) 



After having solved the super-Hamiltonian and super-momentum constraints with no fix- 
ation of the gauge one has to replace 0(i) and <7(i)( r )( s )| r ^ s with their expressions given in 
Eqs.(2.2) and (2.5). The expression in the 3-orthogonal Schwinger time-gauges is obtained 
by putting a (a) (r, a) = 0, y? (a )(r, a) = 0, ? (r, a) = 1 {1) (t, a) = 0. 

With these results we can find the PM expression of 

1) the holonomy along a closed loop T, i.e. P eh A = 
Er=o i> S n>..> Sl >o / ... / A(r( Sl ))..A(r(s n )) d Sl ..ds n , where A[T] = f T A = 
$1 ds A( 7 )( c ) a (rr(s)) r (c) (r (c) are Pauli matrices); 

2) the flux of the electric field across a surface S, i.e. J s d 2 an r 3 e r ^ = E^(S). 
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III. THE PM SPACE-TIME AND ITS INSTANTANEOUS 3-SPACES 

In this Section we illustrate the general properties of PM Einstein space-times and how 
their properties depend on the choice of the York time (selecting a member of our family 
of 3-orthogonal gauges) through the non-local function 3 /C(i) = ■^ 3 K( 1 y After describing 
the proper time of a time-like observer and the properties of the non-Euclidean 3-spaces, we 
compare the results in 3-orthogonal gauges with the IAU conventions in harmonic gauges 
for the Solar System. Then we study the PM time-like geodesies. 



A. The PM Proper Time of a Time-like Observer 

Given a time-like observer located in (r, a r ) (not too near to the particles), the evaluation 
of the observer proper time is done with the line element e ds 2 \( T ^r^ = e 4 g Tr (r, <r r ) dr 2 = 
dl^ T(jr y Therefore from Eqs.(2.14) we get 

d T( T ,&r) = ^e 4 g TT (T, a r ) dr = yjl + 2n (1) (r, a r ) dr = 

\ _ G y- ^y/ggrgyr) K 2 (t) _ 

L C 3 2^ \B-rk(T)\ 1 m 2 c 2 + i^(r) ) 



g> t 3 /C (1) (t,<7)] dr. 



(3.1) 



As a consequence, the proper time depends on the r-derivative of the non-local York time 
3 /C(i) at the position of the observer in the 3-space E T . 



B. The Instantaneous PM 3-Spaces E r 

1. The Spatial 3-Distance on the Instantaneous 3-Space S T 

Let us consider two points on the instantaneous 3-space S T (whose intrinsic 3-curvature 
will be given in Eq.(3.6)) with radar 3-coordinates a r Q and a\. They will be joined by a unique 
3-geodesic f(r, s) = a r Q + (a{-a r ) s + C {1) (r, s), ^(0) = a r Q , f(l) = a[, ^(0) = ^ 1} (1) = 0, 

solution of the geodesic equation S ^ ,s) = - ^2 UV 3 r[ 1)tM) (r, £(r, s)) with the 

3-Christoffel symbol given in Eq.(2.15). 

At order 0(() we get the following solution for the 3-geodesic 



+ « " <) ( f - f ) ds, r ds 2 3 r[ 1)u „(r, a + (a, - B ) s 2 ) 

Jo Jo Jo 



(3.2) 



30 



Since Eqs.(2.14) implies that at the first order the line 3-element joining the two points 

is 



V rs 



ds 



the geodesic 3-distance between the two points is (dE U ciidean(&o,&i) — \&i — <?o\ — 
(°T ~ °~o) 2 1S the Euclidean distance with respect to the flat asymptotic 3-metric) 

d(a , (Ti)(t) = / dS(r) — d EucHdean (a , <7i) + 

J o 

- E W - <) f dsi k (a[ - al) d s (2 + T<?>) - 

- (at - a s ) d r (2 + r?))] (r, <?„ + fa - & ) ■ (3.4) 

As expected it does not depend upon the inertial gauge variable 3 /C(i) 6 . 

Let us remark that in general a 3-geodesic of the 3-metric 3 g(i) rs = — £ A Q(i)rs on the 
3-space S r is not a space-like geodesies of the 4- metric 4 g(i)AB- 



2. The Extrinsic 3- Curvature Tensor 



From Eqs.(1.8) and by using ^aalab = b~ab — §, we get that the extrinsic curvature 
tensor of our 3-spaces in our family of 3-orthogonal gauges is the following first order quantity 

3 K {1)rs (T,a) = <T(l){r){s)\r&>(T,a) + S rs 3 K {1) - d T rW + <9 r n ( i )(r) - ^ <9a™(l)(a)) (t, <?), 

a 

(3.5) 

with fi(i)( r ) and ct(i)(t-)(s) |y^s given in Eqs.(2.4) and (2.5), respectively, and with T-P given 
by Eq.(2.6). Therefore, our (dynamically determined) 3-spaces have a first order deviation 
from Euclidean 3-spaces, embedded in the asymptotically flat space-time, determined by 
both instantaneous inertial matter effects and retarded tidal ones. Moreover the inertial 
gauge variable 3 -^(i) (non existing in Newtonian gravity) is the free numerical function 
labeling the members of the family of 3-orthogonal gauges. The extrinsic curvature tensor 
depends on the local ( 3 K^) and also on the non-local one ( 3 /C(i)) through the shift function. 



6 Instead a space-like 4-geodcsic depends on it. Indeed the extrinsic curvature tensor 3 K rs is a measure, 
at a point in the 3-space S T , of the curvature of a space-time geodesic tangent to the 3-geodesic (3.2) at 
that point, see Refs.[12] 
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3. The Intrinsic 3- Curvature Tensor 



The 3-Riemann tensor is given in Eq.(2.16). The 3-Ricci tensor and the 3-curvature scalar 

are 



3 R(l)rs = 3R W-rus = S rs A (T^ + 2 + d r d s (if) + if > - 2 = 

u 

AttG d r d s {uv) 

- Sr S a (if > + ^ E f r ?>) + * * + r ?> - 5 E I = 

c c 

47rG , 9 r d s (t/y) 1 4 tt 

- ~^3 - V d ™ + -^-J^W(i) + g fy 1 )™' 

3 i? (1) = E 3 %)^ = -8 A0 (1) + 2 E M> = (3.6) 

r a 

We see that, apart from distributional contributions from the particles, the intrinsic 3- 
curvature 3 -R(i) of these non-Euclidean 3-spaces is determined only by the tidal variables, 
i.e. by the PM GW's propagating inside these 3-spaces. 



C. Comparison with the Barycentric Celestial Reference System (BCRS) of 
IAU2000 in the Harmonic Gauge used for the Solar System. 

In Refs.[13] there is the 4-metric chosen in the astronomical conventions IAU2000 to 
describe the Solar System in the Barycentric Celestial Reference System (BCRS) centered 
in its barycenter by using a PN approximation of Einstein's equations in a special system 
of harmonic 4-coordinates x B . The barycenter world-line (a time-like geodesic of the PN 

4-metric 4 gBnv(xB)) is the time axis x b ^(t b ) = (^x° b (t b ); (pj, where Tb is the proper time 

of a standard clock in the solar system barycenter, ((drs) 2 = e gBoo(xB(B)) i.dx° B ) 2 ). It is 
approximately a straight line if we neglect galactic and extra-galactic influences. Through 
each point of this world-line we consider the hyper-surfaces x° B = ct B = const, as instanta- 
neous 3-spaces E x o with rectangular 3-coordinates (practically they are the quasi-Euclidean 
3-spaces of a quasi-inertial frame of Minkowski space-time, even if they do not correspond to 
Einstein's 1/2 clock synchronization convention). In each point of the barycenter world-line 
there is a tetrad with the time-like 4-vector given by the barycenter 4-velocity and with the 
3 mutually orthogonal kinematically non-rotating spatial axes (no systematic rotation with 
respect to certain fixed stars (radio sources) in the instantaneous 3-spaces t B = const.). 
This is a global reference system, with the following PN solution of Einstein's equations 
for the 4-metric 4 (?_b^(xb) (the potentials wb and w B i are static and of order G, so that 
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N B - 3 g%N Bl N B , {x B ) 



gBoo{X B ) 

4 gBoi(x B ) 

A gBij(x B ) = —e 3 g Bi j = —e 



1 - 



2w B 



2w% 



+ 0(cT 5 ) (x B ), 



= -cN Bi (x B ) = -e 



^w Bi 



+ 0(c- 5 ) (x B ), 



[x B 



(3.7) 



Eqs.(3.7) imply an extrinsic curvature tensor 3 K Bi j = j^- (N Bi \j + Nsj\i — d Q 3 g Bi j ) of 
order 0(c -2 ), but the 3-sub-manifolds x° B = const, of space-time (the harmonic 3-spaces) 
are not specified: one has to solve the inverse problem of finding the 3-sub-manifolds with 
the given extrinsic curvature tensor. 

By comparison let us consider the N particles in non-harmonic 3-orthogonal gauges as 
the Sun and the planets of the Solar System. Let us neglect gravitational waves (so that the 
3-spaces have negligible intrinsic 3-curvature except for a distributional singularity at the 
particle locations, [see Eqs.(3.6)], where our approximation breaks down). Then by using 
Eqs.(2.2), (2.3), (2.4), the non-relativistic limit of the 4-metric (2.14) in radar 4-coordinates 
(see the embedding in the Introduction to get world 4-coordinates like the ones of BCRS) 
has the following form 



4 0(i)tt(t, <?) 
A g(i)rr(r, a) 



l-^-^-2^ 3 /C (1) + 0(c- 5 ) (r,a), 



—e 



l g{i)rs{T,o) = -e5 r 



+ ( 9 r 3 /C (1) + 0(c- 5 ))(r,a), 

2 w i 
1 + — + 0(c" 4 ) (r,a), 



w 



W(T 



_ x - 3k-(t) 



G 



W r (T. 



Vi 



V-Vi{T) 



Kir(T) + 



l<?-^( r )| 2 



)■ 



(3.8) 



Also in this 3-orthogonal gauge we can get quasi-static potentials (ignoring the motion 
of the sources and assuming that d T 3 /C(i) and d r 3 /C(i) are slowing varying functions of r) 
and the same pattern as in Eq.(3.7) till the order 1/c 3 included. The main difference is that 
w 7^ w 2 = 0{G 2 ). Here w is the Newton potential and w r the gravito-magnetic one. 

If we choose the special 3-orthogonal gauge 3 /C(!)(r, a) = we recover agreement with 
the Solar System conventions. Let us remark that the instantaneous 3-spaces are not hyper- 
planes due to Eq.(3.5), giving the non-vanishing extrinsic curvature tensor 3 K^ rs = 0(c~ 3 ) 
even if 3 K(i) = 0. 
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See Ref . [14] for the status of knowledge on the possibility of the presence of dark matter or 
of modifications of gravity in the Solar System for explaining effects like the Pioneer anomaly 
(to be mimicked by means of 3 /C(i)(r, a) if needed). Further restrictions on 3 /C(i) near the 
Earth will come from the gravito-magnetic Lense-Thirring (or frame-dragging) effect (see 
Refs.[15], Ref. [16] for Lageos and Ref. [17] for Gravity Probe B) when the experimental errors 
will become acceptable. 

Like in the case of the IAU 4-metric and by assuming that the dependence on the inertial 
gauge variable 3 /C(i) is negligible inside the Solar System, by using the 4-metric (3.8) one 
could reproduce the standard general relativistic effects like the perihelion precession and 
the deflection of light rays by the Sun 7 also in 3-orthogonal gauges. See Ref. [19] for the 
derivation of the Shapiro time delay and for the gravitational redshift induced by the geo- 
potential (by using its multipolar description). With only one body (the Sun) in the limit of 
spherical symmetry one can find the perihelion advance of planets with the standard method 
of using the geodesic equation for test particles (see Refs.[15, 18, 20]). 

D. PM Time-like Geodesies 

Let now us consider a time-like geodesic y^(s) = z fl {a A {s)) = + e A cr A (s) (we use the 
natural adapted embedding of the Introduction) with afnne parameter s and with radar 

4-coordinates cr A (s) = (r(s); o" u (s)j to be used as the trajectory of a planet or of a star. 
The tangent to the geodesic is u^s) = ^f 1 = e^p A (s) with p A (s) = 

At the first order the parametrization of the geodesic (with 4- velocity p A (s)) and the 
geodesic equation are 



° A (s) = ^(s) + ^i)(s) + 0(C 2 ), a A (s) = a A + b A s, 



P ( s ) = = b A + 

as as 



d2 ° ( g ) _ 4p4 ( ( ^ da (s) da c (s) _ 4 A , ( ))b B h c ( o q) 

where cr^(s) = a a + b a s is the flat Minkowski geodesic (with respect to the asymptotic flat 
4-metric). The Christoffel symbols are given in Eq.(2.15). 

The solution of the geodetic equation is 

a A (s) = a A + b A s-b B b c [ dsi [ 1 ds 2 *Tf 1)BC (a + bs 2 ). (3.10) 



^0 



As Cauchy data at s = we take the position y^(0) = x% + e A a A with a A = a A (0) = a A 
and the tangent w M (0) = e A p A (0). 



7 For them a 4-metric approximating the static spherically symmetric Schwartzschild solution is enough: 
see for instance Ref. [18]. 
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For a time-like geodesies the tangent in the origin satisfies eu 2 (0) = 1, i.e. 
e 4 g(i)AB(o"(0))p (0)p B (0) = 1, if the parameter s is the proper time. If -u l (0) = ti l , then we 



have u u (0) = (V 1 + U 2 ;^), U 2 = ^ r (U r ) 2 . Therefore, with b r = U r and with the 4-metric 
of Eq.(2.14), for future-oriented geodesies the condition ew 2 (0) = 1 leads to the following 
result for b A 



V = yJl + U* + d {1) (a ), 
d (i)M = -\l l + U 2 2n ( i)(o- ) - i W r n ( i )(r) (o- ) + 

r 

+ ^(^) 2 (r« + 2 (1) )( ( 7 o )]. 

r 

=► b A = bf o) + 5 A ^d {1) (a ), bf o] = (^/l+U 2 ;U r ). (3.11) 

Therefore, with these Cauchy data and by using Eqs.(2.15), the geodesic and its tangent 
take the form 

t(s) = a T (s) = T +(\Jl+U 2 + d {1) (a )^s- 



- [ dsi / 1 ds 2 ({l+U 2 )d T n {1) + 2\Jl+W y^U u d u n {1) + 
Jo Jo v u 

+ W U " [ - \ ( d u «(!)(«) + dv "(!)(«)) + S uv d T (r« + 2 (1) )] ) (<7 G + U s 2 ) 



d = T(3 K=0) (s) + T(3 K )(s), 



t {3k) (s) = 2 \fl+U*d T 3 !C {1) (* ) -\Ys W r ^ 3 /C (1) ( ( 7 ) - 

r 

- J dsi jf 1 rfs 2 ( - (1 + Z?) <9 T 2 3 /C (1) - 

2 sJl+U 2 UU d u dr 3 /C (1) ~Y WUV d u dv 3 /C (1) ) (a Q + Us 2 ), 
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Si 



cr r (s) = a r (0)+U r s- / ds 1 / ds 2 Ul+W)(d r n (1) +d T n (1){r) ) + 



+ 2\Jl+WY^ UU [ 5 »r d T (r« + 2 - - (9 r n (l){u) - U n (1)(r) )J + 

u 

+ Y, wuv [ 2 tf ™ o v (r« + 2 - <u a r (r« + 2 ) ((t(o) + u s 2 ) = 



def 



O"(3^ =0 )( s ) + °"(3K)( S ) = ^( s Jif=0)(s), 



P A {s) = &f o) +pfi)(s), 



p t ( s ) = vi+w 2 + d ( i)(a )- 

- / ds 2 ((1 + Z?) <9 r n ( i) + 2 \/l + W 2 V W" <9 U n ( i) + 

+ ^ WW [ - i (9 U n (1)(v) + 9 W n (1)(u) ) + <U «9 r (lf> + 2 (1) )] ) (<r + W s 2 ), 

Ml) 

P r (s) =V(r - I ds 2 {{l+U 2 ){d r n {1) + d T n {1)(r) ) + 

+ 2 ^1+W 2 ]T ZT [<5 ur 9 T (r« + 2 - 1 (9 r n (1)(tt) - d u n (1)(r) )] + 

M 

+ Y, UUUV [ 2 5 ™ ^ ( r ^ + 2 ^)) - ^ ^ ( r « ) + 2 ^(i))] ) M°) + w s 2). 

Ml) 

(3-12) 

This is the trajectory of a massive test particle. 

By using Eqs.(2.2) - (2.4), it turns out that all the dependence of the geodesic upon the 

non-local York time is contained in the function T(3^)(s), which contributes with — - to 
the component p T (s) of the tangent. 

Once the time-like geodesic cr A (s) starting at a A = a A (s = 0) and arriving at a A = 
a A (s = 1) is known in terms the 4-metric of Eq.(2.14) and denoted 701, we can evaluate 
the HPM expression of the Synge world function (see Refs. [21-23]), i.e. of the two-point 
function (for a space-like 4-geodesics it has the opposite sign) 
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Q((7 ,(7i) = 



ds ds 



If 1 D da%s) da B (s) 

- dse g AB (a (s)) 

= \f d 

1 J(Toi)0 



l+W{\Jl+W +2p[ 1) {a{s)))- 

- £ W W + 2 P(iM s ^ + (! + (! + 2n (1) (a( S ))) 

r 

- 2 \/l+W 2 ^ U r n (1)(r) - 2 (W r ) 2 (If> + 2 ( i))(<t(s))" 



(3.13) 



This is a 4-scalar in both points (the simplest case of bi-tensors [22]) defined in terms of 
the 4-geodesic distance between them. Its gradients with respect to the end points give the 
vectors tangent to the 4-geodesic at the end points. 
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IV. PM NULL GEODESICS, THE RED-SHIFT, THE GEODESIC DEVIATION 
EQUATION AND THE PM LUMINOSITY DISTANCE 

In this Section we study the null geodesies, the red-shift, the geodesic deviation equation 
and the luminosity distance in PM space-times. 

A. The PM Null Geodesies and the Red-Shift 

Let us now consider a null geodesic y^(s) = z t *(a A (s)) = + t\ <J A (s) through the point 
y v = y M(o) = x % + e^cr A (0) with a A (0) = a A = (r ;a ). It will have the form (4.2) with 
a A = a A . 

However now the tangent vector u^(s) = e^p A (s), with p A (s) = da d ^ = b A — 
b B b c J* ds2 4: T A 1 ^ BC (a + b s 2 ) , is a null vector, e 4 g(i) AB (cr(s)) p A (s) p B (s) = 0. Therefore we 

must require the initial condition e 4 (?(i) ab (o" ) p A (0) p B (0) = e ^g^AB^o) b A b B + 0(( 2 ) = 
on b A = (b T ;b r ). 

By using Eq.(2.14) we get that to each given value of b r there are two values of b T 
determined by the following equation 

[1 +2^^)] (F) 2 -2b T J2 b r n {l)(r) (*o) ~[b 2 + 2j2 (&T (r?> + 20 ( i))K)] = 0, 

r r 

¥ = ±VP + c (1)± (a ), 

C(i )± K) = T Vfr[2n {1) (a ) + (b r ? (r« + 2 <f> {1) )(a )} + \ £ ¥ n {1){r) (a ), 

r r 

b A = b A )± + 5 A ^c (1)± (a ), bf o)± = (±VP;¥). (4.1) 

Therefore we get the following form of a future-oriented null geodesic emanating from a A 
with tangent b A = b A ^ + + 5 At C(i)+(<7 ) 

° A (s) = °o + K)+ + c (i)+) s ~ h (o)+ b fo)+ [ d Sl [ ds 2 4 rf 1)BC (a + 6 (o)+ s 2 ), 

Jo Jo 
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r(s) = r + (V& + c {1)+ (a )) s - 

- [ d Sl f 1 rfJ^ r n (1) + 2^VW"9 u ii( 1) + 
Jo Jo v 

d = r (3x=0) (s) + r ( 3 K) (s), 

r ( 3 X) (s) = - jf tfei J 1 rfs 2 ( - b 2 d 2 3 /C (1) - 

- 2 V 7 ^ £ 6" U d T 3 /C (1) - £ 6" b v d u d v 3 /C (1) ) (<r + b {o)+ s 2 ), 



S PSl 



a r (s) = & r + b r s - J ds 1 J ds 2 (b 2 (d r + d T ^(i)(r)) + 

+ 2 VP bU [ 5 -r dr (r« + 2 - I (9 r n (1)(u) - 9 U n (1)(r) ) 

£ & u 6" [2 5 ru 9„ (if > + 2 (1) ) - <u <9 r (if > + 2 (1) )] ) (a(o) + b (o)+ S2 ) 



+ 



d = °\*K=0)( S ) + ^K)( S ) = ^K=0)( S )i 



This is the trajectory of a ray 0/ light. 
The tangent to the null geodesic is 



p A (s) 



p A (0) 



bf )+pUs) = 

bf o)+ + $ At C(i )+ ((T ) - bf o)+ bf 0)+ / ds 2 4 Tf 1)BC {a + b {o)+ s 2 ), 

Jo 

bf 0)+ + 5 AT c (1)+ (a ), 
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p T (s) = VP + c [1)+ (a )- fds 2 (Pd T n (1) + 2V^ VZTcU (1) + 

+ ^ b u b v [ - i (d u n (1)(v) + d v n ( i )(u )) + (U <9 r (r« + 2 ) (<r(0) + b {o)+ s 2 ) 



6 2 + C(i)+((7 ) +P[ 1)( 3^ =0) (s) +P(i)(S A -)(s), 



P(l)pA-) 



+ 2VPJ2b u d u d T 3 /C {1) + ]T 6" b v d u d v 3 /C (1) ) (<r + 6 (0)+ s 2 ), 



p r (s) = 6 r 



ds 2 ( & (<9 r n ( i) + <9 r n ( i) (r ) ) + 



+ 2 ]T b u [s ur d T (r« + 2 - 1 (9 r n (1)(tt) - <9 U n (1)(r) )] + 

^ 6" 6" [2 5 ru 9, (r« + 2 - 5 U „ 9 r (r« + 2 ) (<r(0) + b [o)+ s 2 ) 



+ 



d - br +P(i)(3^= )( s ) 



(4.3) 



with e 4 g {1)A B(<j(s)) P A (s) p B (s) = + 0(C 2 ). 

The point crj 4 = a A (s = 1) satisfies the equation 



(Tl - T"o) 2 - (<?1 - ^2) 2 = 2 V 6 2 C(i)+((T ) 

- 2 6f o)+ 6f o)+ f d Sl I 1 ds 2 
Jo Jo 



P AT V)bc -J2 b ' r 4r (DBcJ {(To + b {o)+ s 2 ), 

r 

(4.4) 



which gives an idea of the first order deviation of the null geodesic from the flat one joining 
the same two points and a A on the Minkowski light-cone e A r]AB (c A — <J A ) (erf — erf) = 
(a 1 — er G ) 2 = (ti — r ) 2 — J2 r ( a i ~ a o) 2 = 0- See Ref. [24] for the use of a similar equations 
in the IAU conventions for the definition of the radial velocity of stars. 

Let us remark that the already introduced Synge world function f2(cr , 0^) of Eq.(4.5) 
vanishes when evaluated along a null 4-geodesics joining the two points: therefore Q(cr , a) = 
is the equation of the null cone at the point a A . If one solves the equation Q(cr , <xi) = 
in Ti, one can find the emission time transfer function for an electromagnetic signal emitted 
at t in a r Q and absorbed in a[ and then study time delays [23] and their dependence upon 
the York time. 



By using the embedding given in the Introduction we get the following expressions for 
the end points and the tangent vector 
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(4.5) 



With the PM null geodesies one can study the light deflection from a massive body and the 
Shapiro time delay (see for instance Ref.[23, 25]): in both cases the main 3 /C(i)-dependence 
comes from the lapse function n^. 



1. The PM Red- Shift 

x l( T o) 



If Vi(0) = } ° is the unit 4-velocity of the object emitting the ray of light at t g 



and v 2 (s) = 



of the observer detecting it at r s = cr T (s), the emitted frequency w(0), 



the absorbed frequency u(s) and the red-shift z(s) (see Ref.[25]) have the following PM 
expressions 



u(0) 

1 



c^(0)« 1(1 (0) = c^(0) e A p A (0) = cv lfl (0) e A (bf o)+ + 5 M c (1)+ (<7 )), 
ck^(s) v 2 ^(s) =cv 2f i{s) e^p A {s), 



u(s) = v 2fl (s)e^p A (s) 

l + z(s) ' u(0) ^(0)e^(0)' 



z{s) = 1 



Vi M (0) Vb 2 + e?b r 



1 + 



vipip) C(i )+ ((7 ) 



V2»(s) (e£ VP + £ b r ) L v ltM (0) (e? VP + e? b r ) 
v 2 ^{s) c ( i)+(<7 ) +P T {1){ 3 K=0) (s) +p[i )( 3 A - ) (s)] +<P(i)(3 X=0 )(s)) 



(4.6) 



This equation allows to find the dependence of the red-shift z(s) upon the non-local York 
time 3 /C((j(s)). 



B. The PM Geodesic Deviation Equation along a PM Null Geodesic and the PM 
Luminosity Distance 

In the inertial frames of Minkowski space-time the flat null geodesies joining to x 2 
with (xi — x 2 ) 2 = is x^{p) = Xi + (x 2 — x±)p, = = x 2 — x±: this implies 
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\x° - x° 2 \ = a/(^i - x 2 ) 2 = dEudidean{~i-,2), where d E uciidean is the Euclidean spatial distance 
between the two points in the instantaneous inertial 3-spaces. 

In curved space-time we have to solve the equation for the null geodesies (see the previous 
Subsection and the Appendix of the first paper in Refs.[19]). However in astrophysics one 
uses the luminosity distance [25] between the emission point on a star and the absorption 
point on the Earth. We have to find the relation of the luminosity distance with the dy- 
namical spatial distance between the star and the Earth in the dynamical instantaneous PM 
3-spaces. 



1. The PM Geodesic Deviation Equation 

As shown in Ref . [25] , to find the luminosity distance between a point (a star) emitting a 
ray of light (eikonal approximation) and a point (the Earth) where the ray of light (propa- 
gating along a null geodesies) is absorbed, we must solve the geodesies deviation equation 
for nearby null geodesies with the same emission point and propagate the resulting deviation 
vector to the absorption point. 

Let the emitting star S have the world-line Hs{ T {ss)) = %o + e A ( s s0 ( a time-like 
geodesic with parameter ss if the star is considered a test particle) with the unit time-like 

4-velocity Vg(r(ss)) = ^a u s( s s) — e A~^~~ ( s s is the proper time). Let s s = 0, with 
1X5(0) = a A , be the proper time of the emission point. 

Let y»(s) = e' I A a A (s), a A (s) = a A + b A s + cr^(s) be the null geodesic (4.2) followed 
by the emitted ray of light, whose tangent vector k^(s) = e A p A (s), p A (s) = b A + p A ^(s) 

(b A = (VP; b r )), is given in Eq.(4.3). 

At the emission point we have the unit time-like vector u A (a ) and the null vector p A (0) = 
b A + S At C(i)(<t ) satisfying e A g (l)AB (a ) u A (a ) u B (a ) = 1 and e 4 g {1)AB (a o )p A (0)p B (0) = 
8 , respectively. To form a (non-orthogonal) frame at a A we must add two space-like vectors 
E S(\)( a o), A = 1,2 satisfying e A g {1)AB (a ) u A (a Q ) E B (X) = e 4 g (1)AB (a ) p A (0) E B (X) = and 
£ 4 9(i)ab(o- ) E s(\) -^5(Ai) = — $AAi (they span a 2-plane orthogonal the star velocity and to 
the tangent to the ray of light at the emission point). 

A set of four vectors satisfying these conditions is ( 4 r] AB b A b B = 0, 4 r/ AB b A E B ^ S ^ = 0, 

e 4 r) AB b A uf o)s = 1, e 4 g (1)AB (cr ) p A (0) u B (a ) = 1 + (c {1) + - n (m) )(a ) = with 
0Js(cr o ) the emission frequency) 



u A (a ) = u A 0)s -5 AT n {1) (a ), uf o)s = (1; 0, 0, 0), 
p A (0) = b A + S AT c {1) (a ), ^ = (1; 0,0,1), 

E S(\)( a o) = Efo)S(X) + Efl)S(\)( a o)i 

With the 4-mctric (2.14), we have c 4 g {1)AB (Af o) + Afa) (B? o) + Bfc) = A[ o) flfo - £ r Ar {o) B[ o) + 
2 A U BU n w + Afa flfo + BJ o) {A\ o) B r {o) + B{ o) ) n (1)(r) - £ r (A\ o) Bfo + B^) 

2E, A[ ^ (lf>+2 4> (1) ). 
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^(o)s(A) - \ U;e (o)s(A) ,e (o)s(A) ,Uj, 

^(l)s(A)K) = r2 (i)w( (J o)e( o)5(A) ;-(rS 1) +2 (1) )(a o )eJ o)5(A) , 
-(r( 1) +20 (1) )(a o )eJ o)5(A) ,o), 

^ e (o)5(A) e (o)S(Ai) = ^AAi, e (o)5(A) = °- 

(4.7) 



Let the absorbing Earth E have the world-line y E (T(s E )) = rc^ + e^cr^(s E ) (a time-like 
geodesic with parameter s B if the Earth is considered a test particle) with the unit time-like 

A ( \ 

4-velocity v e (t(se)) = ^a u e( s e) = 6 a ( s e is the proper time). Let se = si, with 

a A (si) = crj 4 , be the proper time of the absorption point. 

At the absorption point s = s\ we have the unit time- like vector u E {&i) and the null vector 
p A (s 1 ) = b A +pf 1) (s 1 ), withp^j(si) given in Eq.(4.1), satisfying e A g { i )AB {a 1 ) u E {ai) u E (a l ) = 

1 and e 4 5f(i )A B((Xi) p A (s 1 )p B (s 1 ) = 0, respectively. 

To form a (non-orthogonal) frame at o A we must add two space-like vectors F E , x Ja±), 
A = 1,2 satisfying e A g {1)AB (a l ) u^i) F| (A) = e A g {l)AB (cri) p A (s 1 ) F| (A) = and 
e4 fi'(i)AB(o"i) ^e(a) ^B(Ai) = — ^AAi (they span a 2-plane orthogonal the Earth velocity and to 
the tangent to the ray of light at the absorption point). 

A set of four vectors satisfying these conditions is i^r\AB b A b B = 0, a t\ab b A F^e(X) = ^> 
e 4 fi'(i)AB(o"i)p A (si)M|(o-i) = l + (p[ 1) +n(i)-n(i)( 3 ))((7i) = ^pl with u E ((Ti) the absorption 
frequency) 



«^((7i) = uf o)J5 -<f AT n(i)(<7i), = (1; 0,0,0), 

/( Sl ) = b A +p A ){si) ^ 6 A = (1; 0,0,1), 

B(A)(°"l) = -^(o)B(A) ■^'(l)B(A) (^l) ' 



^(o)E(A) — (^0; /(o)E(A)) /(o)£(A)> 0J ) 

^(Wj^l) = (X) ^(DM^l) /(o)S(A); -(r^ + 20(l))(<Jl) / (o)B(A) ' 
^3 

-(if* + 2 (1) )(a 1 ) / ( 2 o)B(A) , - X P?i)(*i) /(oMA))> 

s^3 

X f(o)E(\) /(o)S(Ai) = ^AAi, /(o)E(A) = 0- 

r^3 

(4.8) 
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We can choose e[ o)sw = f{ o)E(x) = g\ o)(x) with £ r=1>2 9(o)(x) 0( o )(Ai) = ^ and = °- 

As shown in Ref.[25] the deviation vector F M (?/(s)) = e^F A (<r(s)), with y^fff,,) = 0, 
along the null geodesic connecting a A and a A has the following properties: 

A) it vanishes at a A \ 

B) its covariant derivative along the tangent to the null geodesic 



DY A (a(s)) B . . 



ds 



'd B Y A (a(s))+'T A c (a(s))Y c (a(s)) 



is 



(4.9) 

orthogonal to the star velocity u A (a Q ) and to the tangent p a (0) to the ray of light at the 
emission point a A \ 

C) its covariant differential along the tangent to the null geodesic is also orthogonal to the 
Earth velocity u b {cj\) and to the tangent p A (si) to the ray of light at the absorption point 

-A 



or. 



Therefore we have 



Y A {a Q ) = 0, 



DY A (a(s)) 
ds 

D Y A (a( s)) 
ds 



^2 A W E S(\)( a o), 



A=l,2 



B ( A ) f E(A)(' T l)' 

A=l,2 



(4.10) 



The deviation vector is solution of the geodesic deviation equation 



D 2 Y A (a(s)) 
ds~ 2 



= P B (s) (d B [p c {s) (d c Y A {a{s))+*T A D {a{s))Y D {a{s))) 
+ 4 r A BE («(s)) [ P C (s) (d c Y E (a(s))+^ D (a(s))Y D (a(s)))]) 



= 4 g AB (a(s))*R B c D E(*(s)) P c Xs)p u (s)Y»(a(s)) 



(4.11) 



with the initial data Y A (a Q ) = and 



DY A (a(s)) 
ds 



«o = Ea=1,2 A W E S(X)( a o)- 



Its solution, evaluated at the absorption point a A , can be put in the form [25] 
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J2 F E(x 1 )^i)JxAE,S)E s{x)B (a ) 




with 



J A B (E,S) = £ Fi^WJ^EswaM, 



e 4 g il)A c(a 1 )Y A (a 1 )FE iXl) (a 1 ) = - £ J XlX (E, S) 



x 



uJs(cr )/c 



(4.12) 



where us{c ) is the emission circular frequency of the light-ray. The Jacobi map J^ V (E, S) 
maps vectors at S into vectors at E 

2. The PM Luminosity Distance 
The luminosity distance is [25] 

d lum (S, E) = (1 + z) y/\deJj\ = ^4 y/\de~tj~\, (4.13) 

where z is the red-shift of the source as seen by the observer: 1 + z = ws{o~ )/u)e{pi), 
with u>e(cti) the absorption frequency. The corrected luminosity distance is Di um (S,E) = 



In the inertial frames of Minkowski space-time one gets d[ um (S,E) = (1 + 
z) dEudidean(S, E), namely the corrected luminosity distance is the Euclidean spatial dis- 
tance. 



y/\detj\. 



In the weak field approximation, by using cr A (s) = o"( D )(s) + °"(i)( s ) we get 
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Y A (a(s)) = Y A (a {o) (s) + a {1) (s)) = 

= Y A (a (o) (s)) + dYA ^f S)) ^) + 0(C 2 ) = 

= y { A Mo)(s))+Y A (a io) ( S )) + dY{o ff s) K f 1) (s)+o(e), 

0(X (o) 

ObY (*(,))= d< + d< + 

a a y-Ms)) d2Y A^ s)) + g^^j . 
9 c 9 B y ((7(s)) = + + 

By using p A (s) = 6 A as implied by Eq.(4.2), Eq.(4.9) becomes 



DY A (*(8)) _ bB , x , , 5 ,, 



&^)(<T (0 )( S ))+9B^)(<7 (0 )(a)) + 



+ ^^(^(s))^ + ^ 1)bc (<j (o) (s))Y { %(<j {o) (s)) 
+ pfi)(s)d B Y A (a {o) (s)) + 0(C 2 )^ 



+ 



DY A (a {o) (s)) | DY A (a {o) (s)) 
ds ds 

As a consequence, the geodesic deviation equation (4.11) becomes 

Biota a v-<4 . /„\\ i n n v-A . . /„\\ , a a a \sA/ /\\R 



b B b c \d B dcY^ o) (s)) + d B d c Y^(7 io) (s)) + d B dcd + 
+5B 4 rf 1)Cf .(ff (o) ( S )) Y ( %(a (o) (s)) + 2 4 rf 1)BS ( ( r (o) ( S )) ^r^5((7 (0) ( a ))) + 

+^pfl)(s)^9 C ^ ) ( ( T (o) ( S )) = 

= V D4 R(i)DBC E (cT {o) (s))b B b c Y^(a {o) (s)), (4.16) 

with the Christoffel symbols and the Riemann tensor of Eqs. (2.15) and (2.16). 

Therefore we have to solve the following two equations (the dependence upon cr^(s) is 
eliminated by the first equation) 
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b B b c d B dcY^(a {o) (s)) = 0, 



b B b c d B d c Y^{a {o) {s)) = b B b> 



B uC 



4 r] AD ] 



R(i)dbce(?(o)(s)) — 

- d B 4 T( l)CE (a (o) (s))] Y B (a (o) (s))- 

- 2*Tf 1)BE (a (o) (s))d c Y B (a (o) (s))) - 

- b B pf 1) (s)d B d c Y A (a {o) (s)). 



(4.17) 



From Eqs.(4.10) the initial conditions are 



Y&M = Y(1)M = o, 



A=l,2 

D1,:|l,((T, " )( ' s) ^i s=0 = (^^F ( f ) (a (o) ( S )))| s=0 + c (1) K)9 ro F ( 5K)= J] Aa)^)5(a, 



A=l,2 

(4.18) 



Since at the zero order we have Af q AB b A b B = 0, 4 i]a B b A E B ^ S ^ = and e 4 i]A B u f )s ^ B = 
1, due to Eqs.(4.7), the solution of the first equation, satisfying the initial conditions (4.18), 
is 



Y(oMo)(s)) = [e 4 VBcuf o)s (<J? o) (s)-<jC)) J2 A wK)sW 
= (r (0) (s)-r ) A W E t>) 



X=l,2 
9 )Sf(A)> 



A=l,2 



(0) ~ds^~~ = hB 9b Y (°)( a (°)( s ^ = ^2 A m E ?o)s(\)> independently from s. 

A=l,2 

(4.19) 



Let us remark that YA{a{) is proportional to T\ — r — yj {&i — a ) 2 = d EuciideanO-, 0) as 
expected at the zero order in Minkowski space-time. 

Then the second of equations (4.17) and its initial conditions (4.18) become 
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b B b c d B d c YA(a {o) (s)) = {e^ uv uf o)s (al o) (s)-^))b B b c [\ AD a R {1)DBC e{<J {o ){s)) 

~ dB*rfi)CE(<T(o)(s))\ ^2 A W E (o)S(\), 



A=l,2 



(b B B Y&(<7 W (a)))|.=o = A W { E (i)sw(°o) ~ C {1) (<To) E, 



(o)S(A) ) 



(4.20) 



A=l,2 



with E A )s{x) (a ) given in Eq.(4.7). 

Since we have a (o) (s) = a A + b A s, we get b B dY m^ {s)) = Ay^a^s)) and 



■ dY^(a (o) (s)) 
l s=0 ~~ ds l s=0 



Therefore the solution of Eq.(4.17) with the given initial data is 



A=l,2 

+ ds 1 ds 2 
Jo Jo 



Y(iMo)(s)) = [J2 A W (e A )SW (*o) ~ c {1) (a ) Ef o)s{x ^ 

(e 4 VB C uf o)s (<7 ( C 0) (s 2 )-<J C )) 



s + 



b B b c ^ AD 'R {1)DBCE {a {o) {s 2 )) -d B *r( 1)CE (a {o) (s 2 ))) 
Yl A W E (o)s(\)- 



(4.21) 



A=l,2 



By using Eqs. (4.14), (4.19) and (4.20) the last line of Eq.(4.12) becomes 9 



9 We also use b A = (1;0,0,1), e 4 5(1)FA((7l) Ff (Al) (tn) ^ 1)s(A) (a„) = - £ r=1>2 [1 + (if } + 

2 0(l))(ffl)]ff(' o)(Al) 3(' o)(A) , e^lJFA^O-F^jBCAO^IjSCA)^") = £r=l,2 (^ + 2 ^(l))K).9[ o) (A 1 )5('o)(A) 

and cj s (cr )/c = 1 + (c (1) + «(i) - n (1 )( 3 ))(<7 ) = 1 - + + 2 (1) + \ ri(i)( 3 ))(<7 ) (we used Eq.(4.1) 



for c ( i)(<7 )). 
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^ffd^i) y V)4)( ff i) = 



A 



(A) 



J(E,S) Xl \ = ({n -T ) [1 +2 0(i)((7i) + (C(i) - W(i)( 3 ))(<r )] + 

+T(l)((7i) - (C(i) + 2 ( i))((7 o ) Si) 5 Al A + 



)(A) 



r=l,2 

/■SI /-S2 

/ ds 2 / ds 3 (t (o) (s 3 ) - r ) W(i) Ai a(^( )(s3)) = J(E, S) {o)XlX + S)(i) AlA , 

Jo Jo 



3 {$1 E)(o)X 1 X — (Ti — T ) 5 AlA — d Eu clidean(S, E) 8 XlX , 

J(S,E) {1)Xl x = (r ( i ) ( ( j 1 )-(2n (1) +rf ) -in (1)( 3 ) )( ( T )s 1 + 

+ d E uclidean(S, E) (ri(i) + + 2 0(1) - i ™(1)(3))) ^Ai A + 

+ E (dEuclideaniS, E) ^(fTi) - r\ (1) ((Xo) Si) #( o) ( Al ) 0( o 



)(A) 



r=l,2 



/•SI /-S2 

/ C?S 2 / ds 3 (r( )(s 3 ) -T )W(i) AlA ((7( )(s 3 )), 

Jo Jo 



W ( i) AlA ( ( 7( )(s 3 )) = eVi^i)^.) 6 ^" [V D4 %)^cx(^(o)(s3)) - 

-^B 4r fl)CK(°"(o)( s 3)) Ef o)s{x) = 

#(o)(Ai) 6 (^ 4 -R(l)r3rs ~ 4 -R(l)rrrs + 4 -R(l)r33s _ 4 -R(l)rr3s) + 



r,s=l,2 



+ 



(Sr + ft) ( 



4rr I 4-pr 

1 (l)rs + 1 (l)3s 



) 



9( )(\y 



(4.22) 



By using n - r G = d Euc i idean (S, E) we get J( ) XlX (S, E) = d Eu ciidean(S, E) 5 XlX . As a 
consequence we get the following expression of the corrected luminosity distance 
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D lum (S,E) = d ^ {S f ,B) = y/\detJ(S,E)\ = 

1 + Z{Si) 

= \J J (S, E)( )n J (S, -E)( )22 + J (S, E)( )n J (S, -E)(l)22 + J (S, -E)(1)H J (S, -E)( D )22 

, (g f\Ji i J ( S > E h)n + J(S,E) {1)22 n(rA _ 

— dEuclidean{b, £j) \ \ 1 H j (Q ■ h U\^) — 

y U>Euclidean{>5, & ) 



''Euclidean 



(S,E) 



A=l,2 



Ea=1,2 ( ^ ^)(1)AA(3^=0) + J" (5, #)(l)AA(3*0 ) 

\ 1 + d r^E) + 0(C2)? 

^ UEuclidean\>->, & ) 

£ J(S,E) (1)XX = 2(r (1) ( ( 7 1 )-(2n (1) +r( 1) -in (1)(3 ))(a ) Sl + 

+ d E uclidean(S,E) (n(i) + + 2 0(i) - ^ W(i)( 3 ) ) ((7i)) + 

+ (dEuclidean(S,E)ri 1 \a 1 )-Ti 1 \a ) Sl ) (gU^- 

X,r=l,2 

- ds 2 ds 3 (r (o )(s 3 ) - t ) Y W(i)aa(o"(o)(s3)), 
Jo Jo x=l2 

J(S,E) (1)XX(3K) = Sl ((4d T + d 3 ) 3 /C (1) )(<7 ) -irf Mdean (5,E) ((2d T + ^) 3 ^(i))(<7i)- 



A=l,2 



/ ds 2 / ds 3 (r (o) (s 3 ) -r c ) ^ g\ o)[x) g s (o){x) (d r d s 3 IC {1) ) (<r (o) (s 3 )), 

^ - 70 A,r, S =l,2 V ' 



(4.23) 



where Eqs.(2.15) and (2.16) have been used to find the dependence upon the non-local York 



time 3 /Cm- 



Let us remark that Eq. (4.6) implies that the frequency w(0) of the light emitted from 
the star is u(0) = (1 + z(si))u(si), where u(si) is the frequency absorbed on the Earth. 
Since u>(0) = ct>s M (0) e^p A (0) = v rec (S,E) 10 is also the radial (i.e. along the line of sight) 
recessional velocity of the star, we have that the recessional velocity is proportional to the 
red-shift (i.e. it is a red-shift- velocity cz). On the other hand, for small deviations from the 
Euclidean distance, Eq.(4.23) can be written as 

D lum (S, E) » d Euclidean (S, E) + ^J2 J ( S ' E h)>* = « + /?(! + z( Sl )), (4.24) 

A=l,2 



10 



Due to the use of proper time cvg has the dimension of an ordinary velocity with respect to t = t/c. 
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because the term — (2n(i) + I3 — \ n(i)(3))(<r ) si contains w(0), i.e. a linear dependence on 
the red-shift. 

These two results imply that the recessional velocity of the star is proportional to its 
luminosity distance from the Earth (V rec (S, E) = Az(si) + B) at least for small distances. 
This is in accord with the Hubble old redshift-distance relation which is formalized in the 
Hubble law (velocity-distance relation) when the standard cosmological model is used (see 
for instance Ref.[26] on these topics). Again these results have a dependence on the trace of 
the non-local York time, which could play a role in giving a different interpretation of the 
data from super- novae, which are used as a support for dark energy [27]. 
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V. THE PM EQUATIONS OF MOTION OF THE PARTICLES AND THEIR PN 
EXPANSION 

In this Section we study the HPM equations of motion for the particles. We formulate 
the problem of the definition of the center of mass and relative variables first in general 
relativity and then in PM space-times, using the two-body problem as an example. Then 
we study the PN expansion of the equations of motion and we consider the 1PN limit of 
PM binaries for vanishing York time. 



A. The PM Equations of Motion for the Particles 

From Eqs.(5.2) and (5.3) of paper II, by using Eqs.(2.1), we get the following expression 
for the momenta and the equations of motion of the particles (t]^(t), ?)[(t) = 0(1), rrii = 
MO(C), with M the ultraviolet cutoff) 



Ki r [Tj 



r) r (r) M 



dr 



(i + 2(r« + 2 (1) ) 



Hi) ~ E c Viir) [n (1)(c) + (ri 1} + 2 (1) ) fj^r)] 



+ 



n 



(i)W 



l-C(r) 



) 



+ 



dn (1)(a) 9(ri 1) + 2 (1) ) . a 



drjl 



+ 



drfi 



On 



(i) 



(5.1) 



The last line is a consequence of the ultraviolet cutoff, which allows the definition of the 
HPM linearization. 

Eqs(5.1), being implied by Hamilton equations derived from a standard relativistic par- 
ticle Lagrangian (see the action (3.1) in paper I), are equal to the geodesic equations for 
point-like scalar particles notwithstanding these particles are dynamical and not test ob- 
jects (for spinning particles this is not true due to spin- curvature couplings, see for instance 
Ref.[28]). 

Eqs.(5.1) may be rewritten by putting all the terms involving the accelerations at the 



first member. Since Eqs.(2.2)-(2.4) and «j 



rriicr) 



| + M0(() imply that the functions 
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/(i)( r >^) = <l>(i)(T,v),n m (T,a),n( 1){r) (T,a) depend on rf k {r) and rj r k (r) with k = 1,..,N, for 
each of these functions we have £ (r, a) = £^ E s + ^(r) ) . 

Due to the result #(r)=O(0, we get £ / (1) (t,«t) = E s #00 + 0(C 2 ). 

Therefore, the terms involving the accelerations have the following expression 



l-^(r) 



i + 2(r« + 2 (1) ) 



^(D - E c (r) %)( C ) - E c fa c (^)) 2 gg + 2 (1) ) 



#00 



+ 



, „ s/2 [E + 2 ^ r ) ( r ?> + 2 ^))) 

I 1 ~ w)) 3/ c 



+ ro n (i)0) + 

(l-^(r))^ 

+ ^(r) -^(t) (l + 2(r( 1 ) + 2 (1) )-3 

+ -t=v e e w i^r + ™ ( 2 ' (r ^ 20(i)) 



+ 



"(D - E c (r) n (1)(c ) - E c mr) f (lf> + 2 
l-^(r) 



1 — fji(r) u j¥=i 
1 



l-C(r) 



#00 + 



l-C(r) 



9 t); 

g(r c 1} + 2 (1) ) 



)+o(C 2 ). 



+ 



'(r,*(r)) = 



(5.2) 



As a consequence, after having rewritten the lapse and shift functions in the form = 
7i(i) — <9r 3 ^C(i), ^(i)(r) — n (i)(r) + <9 r 3 /C(i), to display their dependence on the inertial gauge 
variable 3 /C(i) = -^ 3 K^, we get the following form of the PM equations of motion of the 
particles 



Vi 



i-^(r) 
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/ _ dn {1) {rMr)) t)[(t) r ^ 



drft 



+ 



9 ™(l)(r) 



l-C(r) 



+ 



u s 

(1) 



dh(i)( u ) 



g(r^ + 2 (1) ) 
0i£ 



- Vi 



7/. r. -L 



(^ c (r)) 2 9(ri 1) +2 (1) ) 



E *M (' + E « ) <*■ «*» - 



€( r ) 
,(1) 



9^ 



+ 



1-Vi(r) 

d ii „ Jr r( r i^( r )iwi( r )) 



+ 

o- 

fl?k 3 ^(D + 2 E #W ^™ + E Vt(r) Vt(r) (r, *(r))) + 0( C 2 ) 



<9ctL 3 }Ct 



d 2 3 /G 



TO,; 



+ o(c 2 ). 



(5.3) 



Since Eqs.(5.3) imply rji ^(r) • ^(r) = rji m^ 1 (1 - r] t (r)) 3/2 ^(r) • ^i(r|rfi(r)|rf^i(r)), the 
final form of the equation of motion of the particles is 



mViVlir) = Vi \l 1 - €{ T ) [fl -V r i( r )Vi(T) ■ ^{r\Vi{r)\VkM T )) = 



def 



ViF r iij\Viij)\Vk^i(r)). 



(5.4) 



The second member of Eqs.(5.4) defines the effective force F[ acting on particle i. It 
contains: 

a) the contribution of the lapse function n^, which generalizes the Newton force; 

b) the contribution of the shift functions which gives the gravito-magnetic effects 

n. 

c) the retarded contribution of GW's, described by the functions of Eqs.(2.6), whose 
contribution at the order 0(C) is given in Eq.(2.7); 



11 Since we have ff^T = ct) = (see Eq.(5.5)), the term EJ a ^- {^g; 



~TTtf THf- J in Eqs.(5.3) is 

proportional to - x Bq, where Bq is the gravito-magnetic field. It is of order (|) 2 as shown in Eq.(5.6). 
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d) the contribution of the volume element 0(i) (0=1 + 6 0(i) + 0(( 2 )), always summed 
to the GW's, giving forces of Newton type; 

e) the contribution of the inertial gauge variable (the non-local York time) 

3 %) = i 3 ^(D- 

Let us remark that in the gravitational case the regularization with Grassmann-valued 
signs of the particle energies leads to equations of motion for the particles of the type 
ViVi( T ) = Vi---- with rjf = with instantaneous act ion- at- a- distance effects coming from the 
lapse ri(i) an d shift n(i)( r ) functions and from the volume 3-element 0(i) = 1 + 6 (f>^y However 
the retardation present in the solution (2.6) for the GW's is not eliminable and formally the 
equations of motion of the particles are of integro-differential type for the N functions t/[(t), 
i = 1,..,N. However, since t)1{t) = 0((), the retardation effects in the GW's are pushed 
to higher HPM order as shown in Eqs.(2.7), so that at the lowest order we have coupled 
differential equations for the particles. This shows that our semi-classical approximation, 
obtained with our Grassmann regularization, of a unspecified "quantum gravity" theory 
does not take into account only a " one-graviton exchange diagram" but also more complex 
structures already present at the tree level (namely they are not radiative corrections) but 
showing up only at higher HPM orders. 

By comparison in electro-magnetism the coupled equations of motion for the charged 
particles and the transverse electro-magnetic field in the radiation gauge, containing the 
Coulomb potential, allow to find the Lienard-Wiechert solution [29] for the transverse vec- 
tor potential with the no-incoming radiation condition. The regularization with Grassmann- 
valued electric charges [7] implies equations of motion of the type i]i(r)=Qi with Qf = 0, 

so that we get the elimination of retardation effects, i.e. Qj7/j(r)(r — \<r\)=Qi Vi( T )- 111 this 
way the difference between retarded and advanced (or symmetric) Lienard-Wiechert solu- 
tions is killed and it is possible to identify the hidden common action-at-a-distance part 
of such solutions and to express the resulting semi-classical Lienard-Wiechert transverse 
electro- magnetic fields in terms of the canonical variables rfi(r), Ki(r) of the particles. This 
implies that the electro-magnetic retardation effects are to be described as radiative correc- 
tions to the one-photon exchange diagram of QED, which is replaced by a Cauchy problem 
with a well defined action-at-a-distance potential. As a consequence, the final equations for 
the particles are second order coupled ordinary differential equations. To reduce the original 
phase space containing the charged particles and the electro-magnetic field in the radiation 
gauge, one added the second class constraints identifying the transverse electro-magnetic 
field with the Lienard-Wiechert solution and one evaluated the Dirac brackets. It turned 
out that the resulting reduced phase space containing only particles has a canonical basis 
spanned by new particle variables r)[, k ir [interpretable as the old ones r/[, K ir , (no more 
canonical with respect to the Dirac brackets) dressed with a Coulomb cloud] with a mutual 
action-at-a-distance interaction governed by the sum of the Coulomb and Darwin potentials. 
In the rest-frame instant form of dynamics [4, 30-32] ,one can find the expression of the inter- 
nal Poincare' generators: p° = Mc, p r ~ 0, j rs , j rr ~ (the potentials appear in the energy 
p° and in the Lorentz boosts J TT ). Then, after having gone from the canonical basis r/[, 

k ir , to a canonical basis containing internal center-of-mass variables = J2i hi 

and relative ones p a , 7r a , a = 1, .., N — 1, (see Eqs. (2.1), (2.2) of Ref.[31]), the rest-frame 

— # 

conditions p r « 0, j TV « 0, eliminated the collective variables: p « 0, if ps f(p a ,^a)- As a 
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consequence, in the reduced phase space there are second order equations of motion only for 
the relative variables p a , ir a . 

In the HPM gravitational case the analogue of the Hamiltonian action-at-a-distance 
Lienard-Wiechert transverse electromagnetic fields are the action-at-a-distance fields 
0(i) (t, a), n(i) (t, a), n (1){r) (r,?), o- ( i )(a)(a) | a ^ 6 (r, ?), of Eqs.(2.2)-(2.5) and the tidal fields 
Ra(r,a) of Eqs.(2.7). To find a reduced phase space containing only particles (how it was 
done in the electro-magnetic case), we have to add the second class constraints which iden- 
tify the tidal fields # s (t, <?) and II a (r, <?) with the HPM solution of Eqs.(2.7) and ((2.9) in 
our family of 3-orthogonal gauges. To get a set of second class constraints for the elim- 
ination of the gravitational field we must add to the existing first class constraints: 1) 
7r^(r, c?) — j^-q 3 -^(i) (t, (?) ps to the super-Hamiltonian constraint written in the form 
0(t, (?) — [1 + 6 0(!)(t, <?)] ~ 0; 2) 9 l (r,a) « to the super-momentum constraints writ- 
ten in the form a (a)(b) | ^>(t, (?) - <7(i)( a )(&) |a^b(r, a) ps 0; 3) n(r, (?) - n ( i ) (r,o ? ) ^ to 
7r n (r, (?) ~ 0; 4) n (a) (r, <?) - n(i)( a) (r, a) ps to 7r n(o) (r, (?) ps 0; 5) /^(r, (?) - -R(i)a(r, (?) ps 
and Ha{f, <?) — ri( 1 ) a (r, ?) ps 0. Differently from the electro-magnetic case, the evaluation 
of the Dirac brackets for the reduced phase space containing only particles at the lowest 
HPM order is trivial, because all the linearized solutions are sums of terms proportional to 
Grrii, i = 1,..,N and the gauge variable 3 /C(i)(r, (?) is a numerical function. Therefore in 
the evaluation of the Dirac brackets of the variables r/[(r), K ir (T), all the extra terms added 
to the ordinary Poisson bracket are quadratic in [Gm 8 Grrijjj^i and can be discarded being 
of order 0(( 2 ) due to the ultraviolet cutoff rrii = M 0((). As a consequence the variables 
tjI, Ki r , are also a canonical basis of the Dirac brackets at the lowest order: the analogue of 
the electro-magnetic Coulomb dressing is pushed to higher HPM order. 



B. The Center-of-Mass Problem in General Relativity and in the HPM Lineariza- 
tion. 

As we have seen Eqs.(5.4), with the gravitational waves T^^r, <?) given in Eqs.(2.7), 
are the final equations of motion for the particles in a reduced phase space containing only 
particles described by the canonical variables f]l(r), /%.(r). The forces appearing in Eqs.(5.4) 
are the gravitational analogue of the electro-magnetic mutual interaction produced by the 
Coulomb and Darwin potentials. 

As said in papers I and II, the instantaneous (non-Euclidean) 3-spaces S T are non-inertial 
rest frames of the 3- universe (an isolated system including the gravitational field in the chosen 
family of space-times) due to the rest-frame conditions P ADM ~ 0. This remains true when 
the gravitational field is expressed in terms of the particles by means of the solution (2.7). 
The gauge fixings to these three first class constraints, eliminating the collective 3-variable 

of the 3-universe, are J A dm = ^adm ~ 0- ^ we introduce the definition K r ADM = — 
(- Eadm) Radmi with R\dm being the gravitational analogue of the (neither covariant nor 
canonical) M0ller 3-center of energy, the conditions K ADM ps imply R r ADM ~ 0. Therefore 
the 3-center of energy is put in the origin of the 3-coordinates on S T (so that there is 
only an external decoupled center of mass of the 3-universe which can be built in terms 
of the ADM Poincare' generators and of the embedding of the 3-space into space-time as 
shown in Refs.[4, 30-32]) and is carried by the reference time-like observer using the radar 
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4-coordinates. In the inertial rest frames of special relativity [4, 30-32] 12 this implies that 
the reference observer has to be identified with the covariant non-canonical Fokker-Pryce 
center of inertia of the isolated system. This is true also in the non-inertial rest frames of 
special relativity as shown in Ref.[4]. As a consequence the isolated system can be seen as 
an external (unobservable) decoupled (canonical non-covariant) Newton- Wigner center of 
mass carrying a pole-dipole structure (the relative motion of the components of the isolated 
system inside the 3-space E T ) identified by the rest mass M and the rest spin S r of the 
system. This is also true in the gravitational case with Mc 2 replaced by Eadm and with 
the rest spin S r replaced by J^dm- The reference observer defining the radar 4-coordinates 
should be replaced also in this case by an ADM Fokker-Pryce center of inertia (a non- 
geodesic time-like observer corresponding to the the asymptotic inertial observers existing 
in our class of space-times, whose spatial axes are determined by the fixed stars of star 
catalogues). 

This is our way out from the the problem of the center of mass in general relativity 
and of its world-line, a still open problem in generic space-times as can be seen from Refs. 
[33, 34] (and Ref.[35] for the PN approach). Usually, by means of some supplementary 
condition, the center of mass is associated to the monopole of a multipolar expansion of the 
energy- momentum of a small body (see Ref.[36] for the special relativistic case). 

Another open problem in general relativity is the replacement of the particle canonical 
variables 7)1 (t), K ir {r) (i = 1, .., N) with relative canonical variables after the elimination of 
the internal center of mass. In the special relativistic electro-magnetic case (see Eqs.(2.1) of 
Ref.[31]) one replaces them with a naif (Newton mechanics oriented) canonical basis t) t + {t), 
K +r (r), p r a {r), 7r ar (r) (a = 1,..,N — 1) with the relative variables in the rest frame being 
Wigner spin-1 3- vectors. The rest-frame conditions imply K +r (r) ~ and their gauge fixings 

def 

K r = J - McR\ « allow to express rf + (r) as a function only of the relative variables, 
rf + rs f r (p a ,n a ), as shown in Refs. [4, 30-32]. In these papers it is also shown how to 
reconstruct the world-lines of the particles by using the embedding of the 3-space S T into 
the space-time. 

However relative variables do not exist in the non-Euclidean 3-spaces of curved space- 
times, where flat objects like fij(r) = ffi(r) — t^(t) have to be replaced with a quantity 
proportional to the tangent vector to the space-like 3-geodesics joining the two particles in 
the non-Euclidean 3-space S T (see Ref. [37] for an implementation of this idea). This prob- 
lem is another reason why extended objects tend to be replaced with point-like multipoles, 
which, however, do not span a canonical basis of phase space. 

However, at the HPM order in our family of space-times we can rely on special relativity 
in the inertial rest frame by using the asymptotic inertial frame (with the asymptotic back- 
ground Minkowski metric) like it has been done in the solution of the constraints and of the 
wave equation for the tidal variables 13 . This allows to define HPM collective and relative 
canonical variables for the particles. Now the collective variables are eliminated with Eqs. 

def 

(2.11) and (2.13), which at the lowest order become pfo and = — M^cR 1 " « like 
in special relativity. If we take into account also the terms of order 0(( 2 ) in Eqs. (2. 11) and 



In these papers there is a complete clarification of the center-of-mass problem in special relativity. 

In the solutions |%(r) — ?7j(r)| is the Euclidean 3-distance between the two particles, which differs by 

quantities of order 0(() from the real non-Euclidean 3-distance on E T as shown in Eq.(3.2). 
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(2.13), we can find the first HPM deviation (depending also on tidal terms) from the special 
relativistic solution and the HPM equations of motion for the relative variables implied by 
Eqs.(5.4). 

As an example let us consider the two-body case in presence of GW's, namely without 
using the retarded solution (2.7). Instead of defining the overall collective variable of the two 
particles and of GW's (like it was done in Ref. [31], where the transverse electro-magnetic 
field was replacing GW's), let us define the naive (Newton mechanics oriented) collective 
and relative canonical variables only of the particles (with masses m\ and m 2 ; M = mi + m 2 
and \x = mi ^ 2 are the total and reduced masses respectively; we choose positive energy 
particles so that we can replace the Grassmann variables 771 and 772 with their mean value 
< rji >= +1, after having done the regularization, as said in footnote 17 of paper I) 



- , \ rn 1 fJ 1 (T)+m 2 f]2 _ , . _ ( _ ( 
Vi2{r) = — , Pi 2 (r) = rjiiV - m{V, 

Kiz\t) = «i(r) + K2(t), 7Ti 2 (r) = — , 

Vi(t) = ff 12 (r) + ^pi 2 (r), ff 2 (r) = ff 12 (r) - ^pi 2 (r), 

Zi(t) = ^±k 12 {t) + tt 12 {t), k 2 (t) = ^Ki 2 (r) -7r 12 (r). (5.5) 

All these quantities are of the type a = a( G ) + + a (2) + 0(( 3 ): the coordinates start 
with 0(1) terms, while the momenta start with 0(() terms due to the ultravilet cutoff. 

The rest-frame condition (2.11) implies K(i)i 2 (t) and a certain expression for /5( 2 )i 2 (t) 
in terms of GW's, ff( ) 12 (r), P( )i 2 (r), n^ 12 (r) and of the non-local York time. 

Instead the condition (2.13) determines 77(0)12 (t) in terms of P( )i2( r ) and 7f( 1 ) 12 (r) 



V(o)i2 (r) « -A {o) (t) p (o )i 2 (r), 

^ J ml c 2 + 7?f 1)12 (r) - ^ J ml c 2 + nf 1)12 (r) 

A (o){r) = = = • 

^m 2 c 2 + nf 1)12 (T) + ^m 2 c 2 + n 2 {1)12 (r) 

(5.6) 

It vanishes in the equal mass case. As a consequence K( 2 )i 2 (t) depends only on relative 
variables. 

Then 7/(1)12(7") is determined in terms of GW's, of particle relative variables and of both 

— * 

the local and non-local York times: 77(1)12 (t) f(i)(T)[rel.var.]. As a consequence the 
particle 3-coordinates depend only on relative variables 
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Vi(r) = ff {o)1 (r) + 77(1)1 (r) « 

~ - A°)( r )) P(o)i2(r) + ^P(i)i 2 (r) + /(^(^[reZ.uar.], 

jfcOO = 77( ) 2 (r) +?7(i)2(r) « 

~ ~(^ + Ao)( r ))p(o)i2W - ^P(i)i2(r) +/(i)(r)[reZ.uar.], (5.7) 

and the world-lines can be reconstructed by using the embedding of 3-spaces into space-time: 

At the lowest order the velocities ff^i(r) depend on P( ) 12 (t), P(o)i2( r )) an d on the function 
A( )(t) and its r-derivative. 

The ADM energy and angular momentum of Eqs.(2.10) and (2.12) become 

Eadm ~ c {^Jm\ c 2 + v?2 1)12 (r) + ^jm\ c 2 + tt 2 1)12 (t)) - 

G *? 1)12 (t) /^/ m i c 2 + ^i)i 2 ( r ) ^/ m i c' + ^W^K 
c3 l^ 12 ^l ^mlc 2 + nf 1)12 (r) + y^m 2 . c 2 + 7f 2 1)12 0V 

,1^2 (r) y^^ 2 / ifV.ftfr)) , r^r^r)) x 
2 (1)12 V A ym 2 c 2 + vf 2 1)12 (r) ^/m 2 c 2 + vf ( 2 1)12 (r) ' 

- £ Kdi 2 (^)) 2 ( y r ° 1)(T, ^ (T)) + _^^t)L=) - 

\/ m i c 2 + ff ( 2 1)12 (r) y/m 2 c 2 + 7f ( 2 1)12 (r) ; 

G y/ m l<? + *ti)i2( T ) V / ^c 2 + 7f 2 1)12 (r) G 7 jr. 
?£^)W')(2j^| + 

1 P(o)12( r )P(o)12( r ) \ 

2 |p (o)12 (r)|3 J + 

+ i^E y rfV[9 r i?,M s - b 9 T i?- b + E^^M s ^ a i?- 6 ])(T,a) + M C 0(C 3 ), 
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We have the following reconstruction of the particle world-lines in the preferred adapted world 4- 
coordinate system defined in the Introduction (f)[(t) = t]\(t = ct)) 



xt (r) = z»(r, ff(r)) - 5?(t) = < + r + e{f r?[(r) = a£ + e£ ct + 
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■ J ADM ~ 



P(o)l 2 ( r ) +P(1)1 2 (^)) ^(1)12^) - (P(o)l 2 ( r ) +P(1)1 2 (t)) 7T(1)12(t) + 

+ P(o)i 2 ( r ) (t( 2 )i2(t) - A°)( T ) K ( 2 )i 2 ( r )) - P(o)l 2 (^) (2)12(1") - A (o) (r) K( 2 )i 2 (t)) 

3 /C ( i)(r,77 (o)1 (r)) 



+ 



+ ^m\c 2 + vf ( 2 1)12 (r) - A (o) (r)) (p\ o)l2 {r) — 
^Ic 2 + ^ 1)12 (r)(^ + A (o) (r)) (p[ o)12 (r) ^- 

/(o 
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)2 " '/(o)2 
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Jd 3 a[j2 W r d s - v s d r ) R- a M- al d T R- b + 2 3 K (1) d r d s (lf> - lf>) + 

+ 2 (9 T if > + d T rf) - ^ E § ^ ^ " ( r > ^ + Mc L °^ 3 )- 

c 



(5.8) 



The dependence of the rest-frame energy on the York time is pushed to order 0(( 3 ). 
Instead the rest-frame angular momentum (or spin) depends on both local and non-local 
York time. 

Let us now consider the equations of motion (5.4). By using Eqs.(5.5) and the previous 
results, they can be written in the following form 



Mr,{ 2 {r) k m[— [A {o) (r)pl o)12 (T) + f {1) (r)[rel.var.]) 

= mr\V(o)i(r)\ff {o)2 (r)) + F;(r\fJ {o)2 (r)\fJ {o)1 (r)), 

(5.9) 
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Since the forces (depending on P( )i2( r ) and ?7(o)j( r )) are of order 0(Q, since A^{r) is 
of order 0(() involving ^^(t) and since P( 1 ) 12 ( r ) is of higher order with respect to P( D ) 12 (' r )) 
Eqs.(5.9) imply the following equations of motion 

HP\ )i2{t) = F{(t\V(o)i(tM(o)2(t)) - ^ F^(r\ff {o)2 (r)\ff {o)1 (r)), 



d 2 

M 



^(l)12( T ) - J~2 ( A (o)( T ) P(o)l 2 (^))_ 

= Fi(T\V(o)i(r)\fj(o)2(r)) +i r 2 r ( r lfe(r)|rf( )i(r)). 



(5.10) 



While the first equation determines the relative motion, the second equation is a consis- 
tency relation connecting the motion of the particle 3-center of mass to GW's and expressing 
the fact that the overall 3-center mass of the isolated system "particles plus GW's" is in free 
motion. 

If the previous discussion is formulated in the reduced phase space containing only par- 
ticles after the elimination of GW's, by forcing them to coincide with the solution (2.7), 
nothing changes except that the second of Eqs.(5.10) should become the time derivative of 
the rest-frame condition 77(1)12 (t) ~ f^(r)[rel.var.]. 

C. The PN Expansion at all Orders in the Slow Motion Limit. 

Due to our ultraviolet cutoff M we have been able to obtain a HPM linearization without 
never making PN expansions. However, if all the particles are contained in a compact set 

of radius l c , we can add the slow motion condition in the form y/e = - c « yj -^p-, i — 1, ..N 

(R mi = 2G J ai is the gravitational radius of particle i) with l c > R M and A >> l c (see the 
Introduction). In this case we can do the PN expansion of Eqs.(5.4). 

Since we have r = ct, we make the following change of notation 

Vi{V = Vi(t), Vi(t) = , (Liit) = — , 

Ur) = «(r) = *g>. (5.11) 

For the non-local York time we use the notation 3 /C(i) (£,<?) = 3 /C(i)(r, <x) = 

By using (1 - x)' 1 ' 2 = J^Zo H* ^ ( valid for ^ < X )> Eqs.(2.7) and (2.2)-(2.4) 

can be written in the form (here = 0(e fc/2 = (^) fc ) is of order | PN) 
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c c 



S (2^-2)!![2(A; 
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Vi<t>(i)(T,ffi(T)) = Vi E ^ ( or m _ - mi + ^ ^(i)i[(2fc)](*>?i(*)l%(*))), 

j# 2 WiW ^jwl fc=i 



(2fc)!! 2|^(0-^(0I 
16 |^(0 -a\ ' ' 



Vin(i){T,f)i(T)) = rji - drlfk 3 /C ( i) (r,rfi(r)) + n (1) (r, ^(r))] 



G m- °° 

+ E ^ -^r 1 E "(i)ii(2fc)](*.?i(*)i%(*)) 

jf^i fe = l 
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Kit) -Vjit) I 
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Vin(i)(r)(r,ffi(T)) = rn d T 3 /C ( i) (r, rf^r)) + n (1)(r) (r, ?^(r) 

= 7}i [<9 r 3 £ (1) (t,?7 i (t)) + Vj ^^«(l)(r)j[(l)](*,?i(*)l%(*) + 

+ ^'-^ (^(l)Wj[(2fc)]+«(l)(r)j[(2fc+l)](*,?i(*)l%(*)) 



fe=l 



"(l)(r)j[(l)](*,?i(*)l%(*)) 



+ 



i(#)-^))f-S(t)-C-W)^ 



2 I^W-i-WI 2 

n(l)(r)j[(2fc+l)] {t,vM\Vj{t)) 



■y 

(2k- 1)!! 



( 2A; ) !! 1^(0 V2 



+ 



•(£(*) -*?,•(*)) ^ ^ft)^ 

2 " '/;(/)-'/,(/) 2 " 



W(l)(r)j[(2fc)](t,?i(<)|%(0) = 



C / 

d 3 ^ 9 lr^tk 



**\tii(t)-vi\ 



d 3 a 2 



871 Iff 1 — (X2 



(5.12) 



All the quantities are even in ^ except the shift functions which have both odd and even 
terms. As a consequence, 4 g(i) TT and 4 g(i) rs are even in ^, but this is not true for 4 (?(i) rr - 

By using Eqs.(5.12), Eqs.(5.4) can be written in the following form after having being 
multiplied by c 2 and rrii (we use (1 — x)^ 1 = YlhLo x> N v& lid f° r x < 1) 
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rmViam = V iFi(t\Vi(t)\Vj(t)) = 



Vimc 



dfl W , n v i Vi dfl (i) , v\ Vj dn (l) 



dfjl 



+ 2- 



c c 



c c 



+ 



(i)H 



(9?7. r c <9ri M 

M " " J^J 



n ^ c c df\\ 



if dh 



(l)(r) 



(9n( 



+ 



u,v 



+£(!)u-4?-£)<^>> 



+ 



Vj d 



c ^ 



* n2 eI-iI« i) + 2 'A«)- 



c 9% 



1 



c c 



+ 



E v ?(*)<(*)^|) ](*,?*(*)) = 



to,- m 



dfjfdfjv 



■3 1 ~ 



+ 



9 77" 
+ 



+ m »E ^ E (^[(2*)] + ^[(2*+l)]) 
jV« fc=l 



where 
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p r _ dn (1)m] v\ v, dn {1)m] ^ v\ v 3 dn (1)j[m 

^ V c df\\ c drft c dfjj ) 

s 1 

and for k > 1 



p r _ ^(i)j[(2fc)] u[ «i dh {1)mk -2)] v\ Vj dn {m{2k _ 2)] 



+ V ( — ^WMP^j] _ ^ (9n (l)(r)j[(2fc-l)] _ ^Tl(i)( r )j[( 2 fe-l)] \ 

ft dhi 



] Vj Vj Vj 0^(l)(u)j[(2fc-3)] | 

c ^ c c dfj? 

Ill) ' L 



u.v 



+ J2 (7) (^(^S[(2fc-2)] + 2 0(l)i[(2fc-2)]) " 
7 ' ^( f S(2 fc -4)] + 20(l) j [(2 fe -4)]) 



2 7 (7 " ^ + ^ " J=r) (rSU_2)] + 20(Di[(2 fe -2)]) + 



pr v-K dn ( i) Wi[(2fc)] _ t£ (9n (1) ( r)j[ ( 2fc)] _ ^ vj_ dk {1){r)mk)] \ 

ii[(2*+l)] V c ft£ c Qfju 2^ c dfj u )^ 

u 11 11 j^i 'J 



V- ^ V- ^ <9n(l)(«)j[(2fc-2)] , v 

+7^-7 — ^ — (5 - 13) 



u,v 



Let us remark that the force F[ contains at the OPN order the Newton force of Newtonian 
gravity. In Eqs.(5.13) it is possible to see that the terms depending on the inertial gauge 
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variable (the non-local York time) 3 /C(i), absent in the Euclidean 3-spaces of Newton gravity, 
are present at the order 0.5PN. Moreover the force F[ contains both even and odd terms 
at all the orders. See Appendix A for the 1PN expression of the ADM Poincare' generators 
with terms till order 0(( 2 ) included. 

In the standard approach in harmonic gauges the first odd terms start at 2.5PN order: 
they are connected to the breaking of time-reversal invariance due to the choice of the no- 
incoming radiation condition and to the effect of back-reaction in presence of gravitational 
self-force with the associated (either Hadamard or dimensional) regularization (see the re- 
view in Ref. [6]). In our approach the Grassmann regularization eliminates the self- force but 
back-reaction is present due to the constancy of the ADM energy and produces the correct 
energy balance for the emission of GW's as shown in paper II. 

Since we are in a non-harmonic gauge, we use a Grassmann regularization and, moreover, 
we are not introducing ad hoc Lagrangians for the particles, it is not possible to make 
comparisons with the standard results known till 3.5PN order [6] (where also the hereditary 
terms are present: we will need higher orders in the HPM expansion to see these terms, if 
them are present in our non-harmonic gauges). 

D. The HPM Binaries at the 1PN Order 

Let us now consider the 1PN two-body problem, which is relevant for the treatment of 
binary systems 15 as shown in Chapter VI of Refs.[6] based on Ref. [40]. 

Since there is no convincing evidence of dark matter in the Solar System and near the 
galactic plane of the Milky Way [41], we shall ignore the 0.5PN terms containing the non- 
local York time in the study of binary systems of stars in some galaxy. Instead in the next 
Section we will see that the non-local York time can be relevant in the 0.5PN simulation of 
dark matter at the level of mass and rotation curves of a whole galaxy. 

If we ignore the York time, the 1PN equations of motion for the binary implied by 
Eqs.(5.13) and (5.12) are (we consider only positive energy, i.e. rji, r] 2 h>< r\\ >, < r\2 >= +1) 



For binaries one assumes ^ ~ \J "7^ << 1, where l c |f | with f(t) being the relative separation after the 
decoupling of the center of mass. Often one considers the case mi w m-i. See chapter 4 of Ref. [6] for 
a review of the emission of GW's from circular and elliptic Keplerian orbits and of the induced inspiral 
phase implying a secular change in the semi-major axis, in the ellipticity and in the period, during which 
the waveform of GW's increases in amplitude and frequency producing a characteristic chirp. If we would 
add terms of higher PN order from Eqs.(5.7), we would get the analogue in the HPM linearization of 
the standard 3.5PN calculations for the inspiral phase before merging and ring-down (see section 5.6 of 
Ref. [6] and Ref. [38] for a review; see also Ref. [39]). Again the Grassmann regularization gives different 
results for the back-reaction. Instead the PM equations of motion Eqs.(5.3) and (5.4) should be used to 
take under control the relativistic recoils during the inspiral phase. 
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Mt)-UtW 



Mt)-UtW 



r 2(lPN) 



(t). 



The last two lines in each equation correspond to the gravito-magnetic force generated 
by the shift functions. These equations identify the 1PN forces F[^ lpN ^(t). 

Let us now reformulate this two-body problem in the canonical basis of center-of-mass and 
relative variables used in Eqs.(5.5)-(5.10) by restricting us to the lowest order ^(r) = fj^t) = 



V(o)i(t)- Since we have /?(i)j(r) 



mi Vj(t) 



7T(l)12(r) = H 



ui(t) 



V2(t) 



=== with «(i)i(r) + /C(i)2(r) = 0, from Eq.(5.5) we get 

"lffl _/m 2 A, J+\\ VeQ( )12(*) 



'i-(^) 



with ff {o) i(r) 



and 



feM = ^1 = -(^+i (o) (i)) where we introduced the definition v (rel)[o)l2 {t) = 

d ^ { °dt^ ■ 111 wr iti n g V(o)i( T ) we ignored the term —~ c dA( ^ P( o )i2(0 : ^ i s °f higher order 
because it depends on 7T(i)i 2 (t). 

Then at the 1PN order we get 



(5 



7T(l)12(r) = fi 



1 + 



l,Ul(f). 



-tf 2 (t) 



1 + 



1 .tf 2 (0' 



]) 



/X^e (o)l 2 (t) ( 1 + \ [(^ - i(o)W) 3 + + i(o)(t)) 3 



%e/)(o)l 2 (^)\ 2 



(5.15) 



Since the 1PN limit of the function A (o) (t) of Eq.(5.6) is ^ m 2 ~? 2) ( 



/J (mi-m 2 ) / ^rei)(o)12(*) \2 



) 2 (so 



that from Eq.(5.6) we have i] (o)12 (t) ~ + O(^)), we get 7T(i)i2(t) = // %e0(o)i2(0 
n|+™f ^ (re0 wi2(t) ^] as the 1PN limit of the re i ative momen tum. 



1 + 



Then from Eqs.(5.8) the 1PN limit of the ADM energy and angular momentum in the 
rest-frame is (at this order there is no dependence on the York time) 



67 



jrs 

J ADM(1PN) 



Eadm(ipn) = ^rriiC 2 + (i (- v 2 rel) ^ o)12 {t) 
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GM 



M 3 1 ' 



1 1 //q x ^reQ(o)12(*) /X , %eQ(o)12(t) P( )l 2 (t) 

2V + M } C 2 + M l C ' |p (o )i 2 (t)| 



\p(o)l2(t)\ 

(P(o)12(0 V (reO (o)12(0 ~~ ^(o)12(0 ^0(0)12(0) 



2M 3 v c 

while from Eqs.(5.10) the equation of motion for the relative variable is 



(5.16) 



d 



V(rel)(o)12(t) o_ 1 fm 2 ~ r m X ~ \ * _ 
It ~ Ji \~M* 1 ( 1PN )W - Jf*2(lPN){t)){t,Vi(o){t),Vi{t)) ~ 

_ r M P(o)12(t) [-,,/-,, o JM V (rel)(o)12( t ) _ _3/i_ ( %eQ(o)12 (g ' P(o)u(t)) 

\p {o) i2(t)\n \ M ) 2M\ |^(o)i2(0l 

GM / 4 _2/i V 



|P(o)l 2 (t)P 



A - %e0(o)12(0 ' P(o)12(0) 

mH^ 12 ^ \^Mt)\ • 



with the forces F[, 1PN Jt) defined in Eqs.(5.14). 

This is the result (ignoring the 0.5PN contribution of the non-local York time; for it see 
next Section) for the 1PN relative motion of binaries in our HPM linearization in the 3- 
orthogonal gauges, where the energy and angular momentum constants of motion are given 
by the corresponding ADM generators (implying planar motion in the plane orthogonal to 
the rest-frame ADM angular momentum). 

Our 1PN equations (5.16) and (5.17) in the 3-orthogonal gauges coincide with Eqs. (2.5), 
(2.13) and (2.14) of the first paper Damour and Deruelle in Ref. [40] (without G 2 terms since 
they are 0(( 2 )), which are obtained in the family of harmonic gauges starting from an ad 
hoc 1PN Lagrangian for the relative motion of two test particles (first derived by Infeld and 
Plebanski [42]) 16 . These equations are the starting point for studying the post-Keplerian 
parameters of the binaries, which, together with the Roemer, Einstein and Shapiro time 
delays (both near Earth and near the binary) in light propagation, allow to fit the experi- 
mental data from the binaries (see the second paper in Ref. [40] and Chapter VI of Ref. [6]). 
Therefore these results are reproduced also in our 3-orthogonal gauge with 3 i^(i) (r, a) = 0. 



This is also the starting point of the effective one body description of the two-body problem of Rcfs. [43]. 
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VI. DARK MATTER AS A RELATIVISTIC INERTIAL EFFECT DUE TO YORK 
TIME 



In this Section we will see that the non-local York time can be relevant in the simulation 
of dark matter at the level of the mass of clusters of galaxies and of the rotation curves of 
galaxies. 

In the first Subsection we consider the 0.5PN equations of motion. In the next two 
Subsections we will consider two of the main signatures of the existence of dark matter in 
the observed masses of galaxies and clusters of galaxies, namely the virial theorem [44, 45] 
and weak gravitational lensing [25, 45, 46]. Then in the fourth Subsection we will reproduce 
the pattern of rotation curves of spiral galaxies [47] . In a final Subsection we discuss which 
information we obtain on the York time from the dark matter data. 



A. The 0.5 Post-Newtonian Limit of the Equations of Motion for the Particles 

At the order 0.5PN, with the non-local York time 4 ^C(i) (with dimension [ 3 /C(i)] = [I] 
since [ 3 -K"(i)] = [^ _1 ]) taken into account, Eqs.(5.13) for the particles become 



d 2 m o 



dt 2 



G 



d 



m-m\ c dt v " 



(1) + 



+ 2 «?(*) 



dfjf 



+ 



(6.1) 



In these equations we can replace the Grassmann variables with their mean value < rji >= 
1, i = 1,..,N, for positive energy particles. 

Therefore at the order 0.5PN the double rate of change in time of the trace of the extrinsic 
curvature, the arbitrary inertial gauge function parametrizing the family of 3-orthogonal 
gauges, creates PN damping terms with damping coefficients 



li{tM)) = + 2 «?(*) d udtk + £ v?(t)vm 9 u d v ) % 1} (t,i(t))- 

u uv 

(6.2) 

For instance the first term corresponds to a damping when <9 t 2 |^. 3 /C(i)(r, ffi(r)) > 0, but it 
is an anti- damping when d%\~ 3 /C(i)(r, ffi(r)) < 0. Since we have [c 2 d 2 3 K^(t, <?)] |a==^( T ) = 
[A d\ 3 K.(i)(t, a)]|- = » the anti-damping (damping) effect is governed by the acceleration 
of the change in time of the convexity (concavity) of the instantaneous 3-space S T near the 
particle as an embedded 3-manifold of space-time. This is a inertial effect, relevant at small 
accelerations of the particle, not existing in Newton theory where the Euclidean 3-space is 
absolute and absent in all the gauges with 3 K (r, a) = (see for instance Ref. [38] for the 
lowest order of PN harmonic gauges). The other damping terms have similar interpretation 
but with an extra dependence on the velocities. 
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We can rewrite Eq.(6.2) in the following form 

li{tM)) = ^%i)(t,m) + 0(( 2 ), (6.3) 

by using ^ = d t +vf(t) d u and by taking into account the fact that the accelerations dvf(t)/dt 
are of order 0(( 2 ). As a consequence Eq.(6.1) can be written in the form 

dt I \ cdt IA n ) dt J 8% jg h \fj.(t) - fj.(t)\ 

+ O(C) (6.4) 

We see that the term in the non-local York time can be interpreted as the introduction 
of an effective (time-, velocity- and position- dependent) inertial mass term for the kinetic 

energy of each particle: m« i— > rrii (l + ^ ^ 3 K,(i)(t, f)i(t)) j in each instantaneous 3-space. 

Instead in the Newton potential there are the gravitational masses of the particles, equal to 
the inertial ones in the 4-dimensional due to the equivalence principle. Therefore the effect 
is due to a modification of the effective inertial mass in each non-Euclidean 3-space: it is 
the equality of the inertial and gravitational masses of Newtonian gravity to be violated! 



B. Dark Matter in Galaxy Masses from the Virial Theorem 

One of the experimental signatures of the existence of dark matter comes from the use 
of the virial theorem for the determination of the mass of clusters of galaxies [44, 45]. For 
a bound system of N particles of mass m (N equal mass galaxies) at equilibrium, the virial 
theorem relates the average kinetic energy < E kin > in the system to the average potential 
energy < U pot > in the system: < E kin >— — | < U pot > assuming Newton gravity. The 
equilibrium condition is supposed to be more valid for clusters of galaxies rather than for 
galaxies (clusters of stars). For the average kinetic energy of a galaxy in the cluster one 
takes < E kin >^ \m < v 2 >, where < v 2 > is the average of the square of the radial 
velocity of single galaxies with respect to the center of the cluster (measured with Doppler 
shift methods; the velocity distribution is assumed isotropic). The average potential energy 
of the galaxy is assumed of the form < U pot >~ —G 11 ^, where M = Nm is the total 
mass of the cluster and 1Z = aR is a "effective radius" depending on the cluster size R 
(the angular diameter of the cluster and its distance from Earth are needed to find R) and 
on the mass distribution on the cluster (usually a ~ 1/2). Then the virial theorem implies 
M ~ H < v 2 >. It turns out that this mass M of the cluster is usually at least an order of 
magnitude bigger that the baryonic matter of the cluster (spectroscopically determined). 

If we consider the 0.5PN limit of the equations of motion for the particles given in eq.(6.1), 
we have to introduce a correction to final form of the virial theorem due to the extra term 
depending on the non local York time. 

Usually the derivation of virial theorem starts assuming that in a self gravitating system 
at equilibrium we have 
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dt 2 

This implies 



r Yl m * i i 2= °- ( 6 - 5 ) 



E mvUt) + E • ^ = 0. (6.6) 



By using Eqs.(6.1) and (6.2) we get 



In the case rrij = m, the mean square velocity is (v 2 ) = vf, and the mean 

gravitational potential energy for particle ( with 1Z = R/2) has the form (U pot ) = 
~jf Ei>j Vj G ^Z^yj -G'^, with M bar = Nm being the baryonic mass. 

Then, if we define ((rj • iij j(t, fj)) = £V ' Vi(t)J 7j(t, ^(i)), we get the following 

result from Eq.(6.7) 

1 I T/2 / -» \ -» 

Therefore for the measured mass M (the effective inertial mass in 3-space) we have 

M = ^ <v 2 >=M bar + ^-((fi-v ) ) 1 (t,fi)) = 

= f M bar + M DM , (6.9) 

and we see that the non-local York time can give rise to a dark matter contribution M DM = 
M - M bar . 



C. Dark Matter in Galaxy Masses from Weak Gravitational Lensing 

Another experimental signature of dark matter is the determination of the mass of a 
galaxy (or a cluster of galaxies) by means of weak gravitational lensing [25, 45, 46]. Let 
us consider a galaxy (or a cluster of galaxies) of big mass M (typically M « M} 2 n ) behind 
which a distant, bright object (often a galaxy) is located. The light from the distant object 
is bent by the massive one (the lens) and arrives on the Earth deflected from the original 
propagation direction. 

As shown in Ref.[25] we have to evaluate Einstein deflection of light, emitted by a source 
S at distance d s from the observer O on the Earth, generated by the big mass at a distance 
dp from the observer O. The mass M, at distance dps from the source S, is considered as a 
point-like mass generating a 4-metric either of the Schwarzschild type (Schwarzschild lens) 
or of the type of Eq.(3.7) (nearly point-like case). The ray of light is assumed to propagate 
in Minkowski space-time till near M, to be deflected by an angle a by the local gravitational 
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field of M and then to propagate in Minkowski space-time till the observer O. The distances 
ds, do, dos, are evaluated by the observer O at some reference time in some nearly-inertial 
Minkowski frame with nearly Euclidean 3-spaces (in the Euclidean case dps = ds — do)- 
If £ = 9d D is the impact parameter of the ray of light at M and if £ >> R s = (the 
gravitational radius), Einstein's deflection angle is a = ^J^ [25] and the so-called Einstein 
radius (or characteristic angle) is 



s dpds V c 2 duds ^ ^ 

A measurement of the deflection angle and of the three distances allows to get a value 
for the mass M of the lens, which usually turns out to be much larger of its mass inferred 
from the luminosity of the lens. 

For the calculation of the deflection angle considers the propagation of ray 

of light in a stationary 4-metric of the type of Eq.(3.8) (with the assumption that <9 T 3 /C(i) 
and (9 r 3 /C(i) are slowly varying functions of time) and uses a version of Fermat's principle 
(see Sections 3.3 and 4.5 of Ref.[25]). In this description the spatial path a (I) (dl = \a\ 
is the Euclidean arc length in an Euclidean 3-space) is the minimum of the variational 
principle (with fixed end points) 5 n dl = 0, where the effective (position- and direction- 
dependent) effective index of refraction is n = 4 g TT + J2 r 4 g Tr If one studies the Euler- 
Lagrange equations for the variational principle, if one ignores gravito-magnetism and if one 
chooses Vr = e [1 - 2 f - 2 d T 3 /C ( i)] as in Eqs.(3.8) with the choice = and with 

the definition 2 d T 3 /C(i) d — — G ^\^\ ■, then the resulting Einstein deflection angle a turns out 
[24] to be the modulus of the vector 



AG M £ 

a = — with M = M bar . + M DM . (6.11) 



Here £ is the impact parameter if the unperturbed spatial trajectory of the ray is 
parametrized as a(l) = £* + I e (e = is the unit tangent vector of the ray satisfying 
e-£ = 0). 

Therefore also in this case the measured mass M is the sum of a baryonic mass Mb ar and 
of a dark matter mass M DM induced by the non-local York time at the location of the lens. 



D. The 0.5PN Two-Body Problem and the Rotation Curves of Galaxies 

To study Eqs.(6.1) in the two-body case (i=l,2), we have to define center-of-mass and 
relative variables in the gravitational case with non-Euclidean 3-spaces, deviating from the 
Euclidean ones by order 0((), at least at the 0.5PN order. 

With the notation of Eqs.(5.5) for the center-of-mass position we put fj[ 2 {t) = V\ )i2^) + 

fj[ 1)12 (t), where fj {o)12 (t) = mMl+piMl (m = m x + m 2 ; /x = is the reduced mass; 

?7i j = 4f (1 + (— )' i+1 \Jl — 2 ^)) is the non-relativistic center of mass and 77m 12 (*) = 0(() is 
a small non-Euclidean correction. 
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The relative position variable is chosen as p\ 2 {t) = Vi{t) ~ V2^) 17 As a consequence we 
have fj[(t) = fjl o)12 (t) +Vl 1)12 (t) + ^p r 12 (t) and fj r 2 (t) = ^ o)12 (t) + f% 1)13 (t) - ^p\ 12 {t)). 

In Subsection VB also a decomposition p r 12 {t) = pl \ 12 (t) + P[i)\ 2 {t) was adopted. This 
decomposition was motivated by the analogous decomposition of the center-of-mass position, 
solution of the equation J t /dm ~ 0- However the equation for the relative motion (see below) 
do not permit to say if the accelerations d 2 p r ^ 12 (t) / dt 2 and d 2 p r ^ 12 {t) / dt 2 are of different 
order or if they are of the same order. Therefore in this Subsection we do not use this 
decomposition. 

As said after Eqs.(5.15), the vanishing of the ADM Lorentz boosts implies 

V( o) i2(t) «0 + O(^). (6.12) 

It eliminates the non-relativistic 3-center of mass by putting it in the origin of the 3- 
coordinates, fj NR (t) 0. Therefore we have fj[ 2 {t) V(i)i2^) = ^(0- 

Then the sum and the difference of the two Eqs.(6.1) gives the following equations of 
motion for the center of mass position fjfo 12 (t) and for the relative variable p{ 2 (t) (we use 
the notation v r (t) = dp r 12 {t)/dt) 

d2 p\m(t) o p\ 2 (t) l dp[ 2 (t) t 

7 + (tJ 12 (t),v(t)) = ^ 7l (t,^p"' 12 (t), t ;(t)) + ^ 72 (t,-^^ 12 (t), t ;(t)), 



u 

+ (^) 2 E^W^^^) 3 %)^-^^)),. (6.13) 



It should be defined as the tangent to the 3-geodesic of S T joining the two points (see the next Eq.(6.16)), 
which is parallel transported along it. See for instance Ref.[37]. At the orders O(C) and 0.5PN the above 
definition is acceptable. 
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We want estimate the contribution of the York time to the rotation curves. Motivated 
by the 1/c factor in front of this term, we can treat the York time term v/c~f + in the 

relative motion equation as a perturbative term added to usual Kepler problem d = — 
G M rf^p • To realize this explicitly we take a solution with circular trajectory such that 
I Pi2(0 \— R — constant and we make a decomposition of the velocity v(t) = v (t) + vi(t) 
such that v (t) = v h(t) is the keplerian velocity such that v Q = \J~^jf — >i?^oo and where 
v\(t) is the first order perturbative correction such that 

-£i = - v ^n{t)i + {tr P i2{timy (6.14) 

at c 

Then we get 

= y dhn^^ihjuih),^)). (6.15) 
At the first order in the perturbation we get 

v 2 {t)=v 2 a ^l-^n(t). jT d*in(ti)7+(*i,pi 2 (ti),u (ti))^ . (6.16) 

Therefore, after having taken a mean value over a period T (the time dependence of the 
mass of a galaxy is not known) the effective mass of the two-body system is 

M eff = ^f^ = M / d *l"(*l)7+(tl,Pl2(*l),V (tl))^) = 

= M bar + M DM . (6.17) 

Again the effective inertial mass in 3-space is the sum of the baryonic matter M^ ar = M 
plus a dark matter term. 

Let us consider the case mi >> m 2 . Let mi be the visible mass of a galaxy and let m 2 
be the mass of either a star or a gas cloud circulating around the galaxy outside outside 
its visible radius. If the 3-space is Euclidean and the Keplerian orbit circular we have that 
the velocity goes to zero when the distance from the galaxy increases, since drke P l ^c(t) _ 

v n{t) = QM-n{t) -^-r^oo 0. Instead from observations one finds that the velocity tend 

to a constant (till where it can be measured) and this so-called problem of the rotation 
curves of galaxies supports the existence of dark matter haloes around the galaxy (see for 
instance Ref.[47] for a review). Again dark matter is an inertial relativistic effect and the 
experimental dark matter distributions can be used to get informations on the non-local 
York time. 

This explanation of dark matter differs: 
1) from the non- relativistic MOND approach [48] (where one modifies Newton equations); 
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2) from modified gravity theories like the f(R) ones (see for instance Refs.[49]; here one gets 
a modification of the Newton potential); 

3) from postulating the existence of WIMP particles [50] . 

Let us also remark that the 0.5PN effect has origin in the lapse function and not in the 
shift one, as in the gravito-magnetic elimination of dark matter proposed in Ref.[51]. 

E. The Non-Local York Time from Dark Matter Data 

As a consequence of the non-Euclidean nature of 3-space in Einstein space-times there is 
the possibility of describing part (or maybe all) dark matter as a relativistic ineriial effect. 
As we have seen the three main experimental signatures of dark matter can be explained 
in terms of the non-local York time 3 /C(i)(r, a), the inertial gauge variable describing the 
general relativistic remnant of the gauge freedom in clock synchronization. 

The open problem is the determination of the non-local York time from the data. From 
what is known from the Solar System and from inside the Milky Way near the galactic plane, 
it seems that it is negligible near the stars inside a galaxy. On the other hand, it is non 
zero near galaxies and clusters of galaxies of big mass. However only a mean value in time 
of time- and space-derivatives of the non-local York time can be extracted from the data. 
At this stage it seems that the non-local York time is relevant around the galaxies and the 
clusters of galaxies where there are big concentrations of mass and the dark matter haloes 
and that it becomes negligible inside the galaxies where there is a lower concentration of 
mass. Instead there is no indication on its value in the voids existing among the clusters of 
galaxies. 

However to get an experimental determination of the York time 3 -^(i) (r, a) = 
A 3 /C(i) (r, a) we would need to know the non-local York time on all the 3-universe: phe- 
nomenological parametrizations of 3 /C(i) (r, a) will have to be devised to see the implications 
for 3 fT( 1 )(r, a). As said in the Conclusions of paper II, a phenomenological determination of 
the York time would help in trying to get a PM extension of the Celestial reference frame 
(ICRS). This would be the way out from the gauge problem in general relativity: the ob- 
servational conventions for matter would select a reference system of 4-coordinates for PM 
space-times in the associated 3-orthogonal gauge. 
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VII. CONCLUSIONS 



In this paper we ended the study of the PM linearization of ADM tetrad gravity in 
the York canonical basis for asymptotically Minkowskian space-times in the family of non- 
harmonic 3-orthogonal gauges parametrized by the York time 3 K(r,a), the trace of the 
extrinsic curvature of the 3-spaces. This inertial gauge variable, not existing in Newton 
gravity, describes the general relativistic remnant of the freedom in clock synchronization: 
its fixation gives the final identification of the instantaneous 3-spaces, after that their main 
structure has been dynamically determined by the solution of the Hamilton equations re- 
placing Einstein equations. It turns out that at the PM level all the quantities depend on the 
spatially non-local quantity 3 /C = ^ 3 K (the non-local York time) with the only exception 
of the ADM Lorentz generators. 

As matter we consider only N scalar point particles (without the transverse electro- 
magnetic field present in papers I and II) with a Grassmann regularization of the self-energies 
and with a ultraviolet cutoff making possible the HPM linearization and the evaluation of 
the PM solution for the gravitational field. 

We studied in detail all the properties of these PM space-times emphasizing their de- 
pendence on the gauge variable 3 /C = 3 K (the non-local York time): Riemann and Weyl 
tensors, 3-spaces, time-like and null geodesies, red-shift and luminosity distance. Also the 
Ashtekar variables of PM space-times were evaluated in the York canonical basis. We also 
studied the PM equations of motion of the particles, their PN expansion and the PM problem 
of the determination of center-of-mass and relative variables for the particles. 

All the main measurable quantities turn out to have a dependence on the non-local York 
time. However it seems plausible that inside the Solar system and in every nearly isolated 
binary system this gauge quantity is negligible. This is not true at the astrophysical level 
especially for galaxies and clusters of galaxies with a big mass. 

In Section VI we have shown that the main features of the experimental signatures for dark 
matter (masses of clusters of galaxies, rotation curves of spiral galaxies) can be explained in 
terms of the non-local York time at the 0.5PN level in the PN expansion. 

This opens the possibility to explain dark matter inside Einstein theory without modi- 
fications as a relativistic inertial effect: the determination of 3 /C(i) from the mass and the 
rotation curves of galaxies [45, 47] would give information on how to find a PM extension 
of the existing PN Celestial reference frame (ICRS) used as observational convention in the 
4-dimensional description of stars and galaxies. 

Therefore what is called dark matter would be an indicator of the non-Euclidean nature 
of 3-spaces as 3-sub-manifolds of space-time (extrinsic curvature effect), whose internal 3- 
curvature can be very small if it is induced by GWs. It is the Newtonian equality of 
inertial and gravitational masses in Euclidean 3-space to be violated, not their equality in 
the 4-dimensional space-time implied by the equivalence principle. 

This conclusion derives from the analysis of the gauge problem in general relativity done 
in the Conclusions of paper II. The gauge freedom of space-time 4-diffeomorphisms implies 
that a gauge choice is equivalent to the choice of a set of 4-coordinates in the atlas of the 
space-time 4-manifold and that the observables are 4-scalars. At the Hamiltonian level the 
gauge group is deformed and the Hamiltonian observables are the Dirac observables (DO), 
which generically are only 3-scalars of the 3-space. However, for the tidal variables and the 
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electro-magnetic field there is the possibility (under investigation by using the Newman- 
Penrose formalism [52]) that 4-scalar DO's describing them could exist. 

On the other side at the experimental level the description of baryon matter is intrinsically 
coordinate- dependent, namely is connected with the conventions used by physicists, engineers 
and astronomers for the modeling of space-time. As a consequence of the dependence on 
coordinates of the description of matter, our proposal for solving the gauge problem in our 
Hamiltonian framework with non- Euclidean 3-spaces is to choose a gauge (i.e. a 4-coordinate 
system) in non-modified Einstein gravity which is in agreement with the observational con- 
ventions in astronomy. Since ICRS [5] has diagonal 3-metric, our 3-orthogonal gauges are 
a good choice. We are left with the inertial gauge variable 3 /C(i) = ^ 3 -^(i) not existing 
in Newtonian gravity. As already said the suggestion is to try to fix 3 /C(i) in such a way 
to eliminate as much dark matter as possible, by reinterpreting it as a relativistic inertial 
effect induced by the shift from Euclidean 3-spaces to non-Euclidean ones (independently 
from cosmological assumptions). As a consequence, ICRS should be reformulated not as a 
quasi-inertial reference frame in Galilei space-time, but as a reference frame in a PM space- 
time with 3 -ft'(i) (i.e. the clock synchronization convention) deduced from the data connected 
to dark matter. Then automatically BCRS would be its quasi-Minkowskian approximation 
(quasi-inertial reference frame in Minkowski space-time) for the Solar System. This point 
of view could also be useful for the ESA GAIA mission (cartography of the Milky Way) [53] 
and for the possible anomalies inside the Solar System [14] . 

Moreover our approach will require further developments in the following directions: 

a) Find the second order of HPM to see whether in PM space-times there is the emergence 
of hereditary terms (see Refs.[6, 54]) like the ones present in harmonic gauges. Like in 
standard approaches (see the review in Appendix A of paper II) regularization problems 
may arise at the higher orders. 

b) Study the PM equations of motion of the transverse electro-magnetic field trying to 
find Lienard-Wiechert-type solutions (see Subsection VB of paper II). 

c) Dark energy in cosmology [27]. Take a perfect fluid as matter in the first order of HPM 
expansion [55] adapting to tetrad gravity the special relativistic results of Refs.[56]. Since in 
our formalism all the canonical variables in the York canonical basis, except the angles 9 l , 
are 3-scalars, we can complete Buchert's formulation of back-reaction [57] (see also Ref.[58]) 
by taking the spatial average of all the PM Hamilton equations in our non-harmonic 3- 
orthogonal gauges. This will allow to make the transition from the PM space-time 4-metric 
to an inhomogeneous cosmological one (only conformally related to Minkowski space-time at 
spatial infinity) and to reinterpret the dark energy as a non-linear effect of inhomogeneities. 
The role of the York time, now considered as an inertial gauge variable, in the theory of back- 
reaction and in the identification of what is called dark energy 18 is completely unexplored. 

Let us remark that in the Friedmann- Robertson- Walker (FLW) cosmological solution the 
Killing symmetries connected with homogeneity and isotropy imply (r is the cosmic time, 
a(r) the scale factor) 3 K(t) = — = —H, namely the York time is the Hubble constant. 
However at the first order in cosmological perturbations we have 3 K = —H + 3 K^- ) with 



As we have seen the red-shift and the luminosity distance depend upon the York time, and this could play 
a role in the interpretation of the data from super-novae. 
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3 i^(i) being an inertial gauge variable. Instead in the spatial-averaging method of Ref. [57] 
one gets that the spatial average of the York time (a 3-scalar gauge variable) gives the 
effective Hubble constant of that approach. Therefore we will try to see if also dark energy 
can be considered as an inertial effect of the York time 19 in the transition from astrophysics 
to cosmology. 



In paper I we showed that in the York canonical basis the York time contributes with a negative term to 
the kinetic energy in the ADM energy. It would also play a role in study to be done on the reformulation 
of the Landau-Lifschitz energy-momentum pseudo-tensor as the energy-momentum tensor of a viscous 
pseudo-fluid. It could be possible that for certain choices of the York time the resulting effective equation 
of state has negative pressure, realizing in this way a simulation of dark energy. 



78 



Appendix A: The PN expansion of the Weak ADM Poincare' Generators. 

The 1PN expansion of the ADM Poincare' generators of Eqs. (2.10)-(2.13) with 0(C) 
order included is 
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Since we have Ra,T ( a ] = #0(4) the GW kinetic term in the ADM energy is of order 
Mc 2 0(\). The expression of the energy has been obtained by making an integration over 
3-space and by using the integral given after Eq.(4.17) of paper II. 
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